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PERTURBATION SOLUTIONS FOR THE INFLUENCE OF 
FORWARD FLIGHT ON HELICOPTER ROTOR 
FLAPPING STABILITY 
By Wayne Johnson* 

U. S. Army Air Mobility RflD Laboratory 
Moffett Field, Calif,, 9403S 

SUMMARY 


The stability of the flapping motion of a helicopter rotor blade in 
forward flight is investigated, using a perturbation technique which gives 
analytic expressions for the eigenvalues, including the influence of the peri- 
odic aerodynamic forces in forward flight. The perturbation solutions are 
based on small advance ratio y (the ratio of the helicopter forward speed to 
the rotor tip speed). The results are valid to approximately y * 9 5, which 
covers the forward speed range of most helicopters. The rotor configurations 
considered are a single, independent blade; a teetering rotor; a gimballed 
rotor with three, four, and five or more blades; and a rotor with N 
independent blades. The eigenvalues of a constant coefficient approximation 
to the flapping equation are obtained by an expansion in y, and are compared 
with the perturbation solution including the periodic coefficients. The con- 
stant coefficient approximation with the equations and degrees of freedom in 
the nonrotating frame represents the flap dynamic reasonably well for the 
lower frequency modes, although it cannot, of course, be completely correct. 
The transfer function of the rotor flap response to sinusoidal pitch input is 
examined, as an alternative to the eigenvalues as a representation of the 
dynamic characteristics of the flap motion. 


INTRODUCTION 


The fundamental motion of helicopter rotor blade is flapping motion: 
first mode out of plane (vertical) displacement from the plane of rotation of 
the rotor. For an articulated rotor this motion is rigid body rotation of the 
blade about a hinge at or near the center of rotation. For a cantilever rotor 
it is elastic bending motion, primarily about flexibility at the blade root. 
Due to the high centrifugal forces on the blade, the mode shape for first mode 
bending of a blade with cantilever root restraint is nearly the same as for 
the rigid body flap motion of a hinged blade; it is the difference in the 
natural frequencies of the cantilever and articulated blades; that is, of 
primary significance for the flap dynamics. This flapping motion of the rotor 
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has a basic role in helicopter stability and control, blade loads, vibration, 
and in most other areas of helicopter behavior. An important aspect of heli- 
copter rotor flapping dynamics is the dynamic stability of the motion in for- 
ward flight. Besides the question of high speed stability, this topic also is 
of interest for its implications on the decay of transient motions of the 
blade, for example, in response to control inputs or aerodynamic gusts. Any 
significant degradation of the blade stability and response to control due to 
forward speed will have important effects on all aspects of the helicopter 
behavior. 

In hover (zero forward speed) , the aerodynamic forces on the rotor blade 
provide the flap motion with high damping, hence, good stability and fast 
response. The motion of the blade in hover is describable by constant coeffi- 
cient linear differential equations, which may be analyzed to obtain the 
dynamic behavior by the standard methods of linear system theory. In forward 
flight however, the aerodynamic forces on the blade introduce time varying - 
specifically, periodic around the rotor azimuth - coefficients into the linear 
differential equations describing the flap motion. These periodic coeffi- 
cients are due to the once-per-revolution variation of the free stream veloc- 
ity seen by the rotating blade when the helicopter has forward speed. Systems 
described by periodic coefficient equations require a considerably more 
involved analysis in order to investigate their stability and response to con- 
trol. Appendix A discusses the behavior characteristic of periodic coeffi- 
cient systems, and presents some of the results of the mat’ ematical theory of 
such systems. Often a direct numerical integration of the equations of motion 
is used as either the most or the only practical solution method; the mathe- 
matical theory of periodic coefficient equations has also been used for the 
analysis of rotor dynamics in forward flight. 

The flapping stability of a rotor in forward flight has been investigated 
in a number of publications (refs. 1-20). The primary reason for the frequent 
reexamination of this one problem is the search for a satisfactory technique 
for the treatment of the influence of the periodic coefficients in the differ- 
ential equations, for even to determine the stability of such an equation is a 
much more difficult mathematical problem than with constant coefficients. 
Usually the flap stability has been determined by numerically integrating the 
equations of motion, and then using the results of Floquet theory to find the 
eigenvalues (roots) from the transient solution over one period (refs. 2, 4, 

5, 7, 10, 12, 13, and 18; the references usually describe the mathematical 
theory required for the method) . There have also been a number of solutions 
using the methods of perturbation theory (refs. 3, 6, 9, 14-17), using an 
analog computer to solve the equation (refs. 8, 11, and 20), and a solution 
using the classical methods of the analysis of Hill's equation (ref. 1). The 
reverse flow region of the rotor is important in very high speed forward 
flight, and several investigations have considered the aerodynamics of reverse 
flow and its impact on the flapping stability (refs. 5, 7, 8, 11, 12, 14, 15, 
and 19). Pitch/flap coupling of the blade motion is another import factor 
which has been treated (refs. 2, 3, 8, 12, 14, 15, and 17). A number of the 
studies have extended the stability analysis by adding other degrees of 
freedom besides rigid flapping to the blade motion description: flapwise 

elastic bending of the blade (refs. 4, 5), blade torsion (refs. 5, 9, 11), 
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blade lag motion (refs. 10, 16), or the flap motion of the other blades of the 
rotor (ref. 13). There is, however, still work that may be done with just the 
problem of flapping stability, particularly since most of these studies have 
only presented numerical results. 

This report investigates the effect of forward speed on the helicopter 
rotor flapping stability and response to control. A solution for the stability 
is obtained by a perturbation technique which is reasonably straightforward, 
and which provides analytic expressions for the eigenvalues of the flap 
motion. The solution method is direct enough so that it is possible to con- 
sider here a number of rotor configurations, in addition to the problem of a 
single independent blade which is the subject of most of the literature. The 
perturbation method to be used here is known as the method of multiple time 
scales. It is described briefly in Appendix B, but in fact the method is best 
discussed by example, of which there will be several here. For more informa- 
tion on the mathematics of this perturbation technique and others, the reader 
is directed to reference 21 (and to ref. 14, which discusses the techniques 
specifically in the context of the problem of rotor flapping stability). For 
the speed range of most helicopters, the advance ratio y (the ratio of the 
helicopter forward speed to the rotor tip speed) is a small parameter; a maxi- 
mum of y = 0.4 to 0.5 may be assumed. Therefore, a perturbation solution 
based on small y may be expected to be applicable over the entire range of 
interest for most helicopters. Most of the present work will be concerned 
with the small y results, although the stability at high u (around 1.0 to 
2.5) will be briefly discussed. For the range of speed involving small 
advance ratio, y * 0.0 to 0.5 or so, a perturbation solution to order y 2 is 
satisfactory. It is therefore possible to neglect the effects of the reverse 
flow region on the rotor (since they are of higher order than y 2 - see 
ref. 8) , which results in a considerable simplification of the aerodynamic 
forces which must be considered. 

With the great difficulties involved in the mathematical analysis of 
periodic coefficient differential equations, it is natural to consider the use 
of a constant coefficient approximation to the equations describing the sys- 
tem. That is, the periodic coefficients are replaced by their average values, 
so the equations of motion are reduced to constant coefficient equations which 
may be easily analyzed by the standard techniques. The validity and applica- 
bility of the constant coefficient approximation to several of the rotor con- 
figurations considered here will be investigated. An expansion of the 
eigenvalues for small y will be found, for comparison with the perturbation 
solutions including the periodic coefficients. 

The following topics will be considered in this report. First, the 
stability of the flap motion of an independent blade will be investigated; 
that is, the case of a rotor with a fixed shaft, so that the motion of any one 
blade is independent of that of the other blades of the rotor. A perturbation 
solution (for small y) will be found for the influence of forward flight on 
the eigenvalues. This analysis will be followed by a summary, which will 
collect the results for the eigenvalues. The behavior of these results will 
be discussed, and examples given of their application. Next the flap stabil- 
ity of a teetering rotor will be investigated, followed by a summary and 
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discussion of the results. Then a gimballed rotor will be considered, for the 
cases of three, four, and five or more blades. The description and behavior 
in the nonrotating frame of a rotor with N independent blades will be given; 
it is more appropriate to make the constant coefficient approximation in the 
nonrotating frame since some of the influence of the rotating periodic coeffi- 
cients is retained. Next, the constant coefficient approximation in the 
rotating and nonrotating frames will be considered, and the results compared 
with the solution including the periodic coefficient influence. Then the 
transfer function, that is, the response of the flap motion to sinusoidal con- 
trol inputs, will be investigated, including the influence of the periodic 
aerodynamic forces in forward flight. The transfer function is an alternative 
to the eigenvalues as a representation of the dynamic characteristics of the 
system. Finally, some of the previous work with this problem will be 
discussed. 
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NOMENCLATURE 

constant in the order y solution for g 
constants for the 1/rev critical region 

secular equation parameter, defining the root behavior near 
the critical regions 

transfer function 

/Tf 

pitch/flap feedback gain, Kp = tan 63 
blade index, m ■ 1 , . . ., N 
flap moment due to flapping velocity g 
flap moment due to flapping displacement g 
flap moment due to pitch control 0 
number of blades 
blade radial station 
blade radius 

blade flap degree of freedom 

coefficients in expansion of g as series in y 

rotor tip path plane pitch degree of freedom 
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rotor tip path plane roll degree of freedom 
rotor coning mode degree of freedom 
vector of Bic^is degrees of freedom 
lock number of the blade 

coefficients in expansion of y as series in y 

pitch/flap coupling parameter 

denominator of the hover transfer function 

y ' r 0 

air density 

blade flap mode shape 

blade pitch control 

eigenvalue 

order 1 eigenvalue (the hover limit) 

advance ratio, helicopter forward speed divided by rotor tip 
speed 

y at boundary of critical region 

corner y for 1/rev critical region 

blade flap natural frequency (rotating, per rev) 

blade azimuth angle; nondimens ional time variable 

time scales, i|> n ■ y n i|> 

azimuth angle of the mth blade 

frequency 

rotor rotational speed 

complex conjugate; for the transfer functions, magnitude of 
the input and response 

time derivative (when dimensionless, the derivative with 
respect to ♦) 



The equations and parameters used in this analysis are dimensionless, 
based on the air density p, the rotor rotation speed ft, and the rotor 
radius R. 


PERTURBATION SOLUTION FOR HELICOPTER ROTOR FLAPPING STABILITY 


Equation of Motion 

Consider the first mode flap motion of a single blade of a helicopter 
rotor. The rotor shaft is fixed so there is no coupling of the blades through 
the shaft or control system. Then each blade is independent of the motion of 
the other blades. The equation describing the flap motion in the rotating 
frame is 


e + v 2 3 = y[m- 3 + M g 3 + M 0 (0 - K P 3)J 


This is the equation for small perturbations of the flap motion from a trim 
state. Only blade pitch control (also a perturbation from the trim state, 
that is, cyclic and collective control as required for the given thrust and 
forward speed) is included as an input. The degree of freedom representing 
the perturbed flap motion is 3. For an articulated rotor 3 is the angle of 
rotation of the blade about the flap hinge. In general, (i.e., for a canti- 
lever rotor, or an articulated blade with flap hinge offset) the blade out of 
plane deflection is Bn, where n(r) is the mode shape of first mode flapping, 
normalized to unity at the tip (at r = 1) . The blade pitch control perturba- 
tion is 9; it is input by the control system, as rigid pitch motion of the 
blade about a feathering axis at the blade root. So 6 is the pitch change 
of the blade all along its span. The dot denotes the derivative with respect 
to the blade azimuth angle Jj», which is the nondimens ional time variable. 

The first term on the left-hand side of equation (1), 3, is the flapping 
inertia. The second term, v 2 B, is the flap spring, which has structural and 
centrifugal contributions. The parameter v is the rotating natural fre- 
quency of the flap motion (per rev since frequencies are nondimensionalized by 
ft). For an articulated blade, that is, a flap hinge with no offset or spring 
restrain, v = 1; this is entirely the centrifugal stiffening. A blade with 
structural restraint (cantilever root, or a hinge spring), or with the flap 
hinge offset from the center of rotation, has v > 1; v a 1.1 to 1.15 is 
typical of a cantilever rotor, and \> = 1.03 to 1.05 for an articulated rotor 
with hinge offset. The right-hand side of equation (1) is composed of the 
aerodynamic flap moments, due to the flapping velocity and displacement, and 
the blade pitch control. Mechanical pitch/flap coupling is included, tha* is, 
feedback control of the form A0 * -K p B. This feedback is usually accom- 
plished (for articulated rotors) by mechanical arrangement of the flap hinge 
and pitch form equivalent to rotation of the flap hinge by the angle 5 3 , Sv> 

Kp * tan 63. Positive K p (i.e., 63 > 0) is negative feedback, which intro- 
duces a positive aerodynamic spring into the flap equation through the 
coefficient Mq. 


1 

* 
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The factor y is the rotor Lock number, defined as y = pacR 4 /Ib (where 
0 is the air density, a is the blade section lift-curve slope, c is the 
blade chord, R the rotor radius, and 1^ is the moment of inertia of the 
flap motion) . The Lock number represents the ratio of the aerodynamic forces 
to the inertia forces on the blade; y typically has values from 6 to 10 for 
current helicopters. The left-hand side of equation (1) is the blade aero- 
dynamic flap moments due to flapping and pitch control. Neglecting reverse 
flow, and assuming the mode shape is n 3 r, then the aerodynamic moments are 

m b = -(* + Z u sin *) | 


M 


8 


-y cos 



1 

J 


y sin 


*) 


( 2 ) 


M e = -| + J y Sin \|» + J (y sin i |») 2 j 

where u is the rotor advance ratio . the ratio of the helicopter forward 
speed to the rotor tip speed; for current helicopters, the advance ratio at 
maximum forward speed is typically between 0.3 and 0.5. The neglect of 
reverse flow is valid to about y = 0.5, for reverse flow adds order y 4 
terras to the harmonics of the flap moments. Hence, neglecting reverse flow is 
consistent with the present small y perturbation analysis, which will be 
carried to order y 2 . 


The equation of motion then, substituting for the flap moments and 
dropping the pitch control forcing terms to obtain the homogeneous equation, 
becomes 


* 



I 




8 ♦ ♦ -fa y sin ^8 + |v 2 + y cos y sin 

♦ Kpjjj- ♦ £ y sin <J> + £ (y sin ^) 2 JJfi * 0 (3) 

The derivation of this equation may be found in the literature, for example 
in reference 8, which also considers the reverse flow aerodynamics. The only 
parameters are the flap natural frequency v, which is always 1 or slightly 
above; the pitch/ flap coupling Kp, which is frequently zero for helicopter 
main rotors; the Lock number y, a measure of the relative strength of the 
aerodynamic forces on the blade; and the forward speed y. For the hover 
limit, y = 0, equation (3) reduces to a constant coefficient linear differ- 
ential equation. For forward flight, n introduces periodic coefficients into 
the differential equation. The stability of a system is defined by its eigen- 
values, or roots; there are two for this sr.cond-order equation. The stability 
of the flap motion will be examined for ti;e case of small advance ratio, by 
means of a perturbation technique to hand - <i the influence of the periodic 
coefficients due to forward flight. 
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Hover 


In the hover limit, p * 0, the equation reduces to 

B ♦ $ 8 ♦ (v 2 ♦ Kp * 0 (4) 

The eigenvalues X of this constant coefficient differential equation are 
obtained from the characteristic equation 

X 2 + J X + v 2 + K p J = 0 


which has the solution 



These roots are usually complex conjugate pair, that is, the radicand is 
positive. For that case, the two roots lie on a circular arc in the X plane, 
with center at ReX = -Kp and radius ♦ Kp'fc, so the circle goes through 
the imaginary axis at ±iv. The real part of the roots is given by y only, 

ReX ■ -y/16, which then determines the location of the two roots on this arc. 
The root locus for varying y then consists of the portion of this circle to 
the left of the imaginary axis (since y is positive) , plus the real axis 
from -oo to -Kp. At y » 0 the roots are at X * ±iv on the imaginary axis. 
For y -*■ co the roots approach X ■ -oo and X = -Kp (which notice is t he cente r 
of the circle). The y locus intercepts the real axis at X = -Kp - /v 2 + Kp^ 
(the center of the circle plus its radius), when y/16 = Kp + A>2 + Kp2. For 
y still larger there are two real roots, that is, the radicand of equa- 
tion (5) is negative; for smaller y the roots are a complex conjugate pair. 
For Kp = 0 this intercept occurs at y/16 ■ v, hence at y approximately 
16, which is quite large for current helicopters. Negative pitch/ flap 
coupling Kp (positive feedback) is required for the intercept to occur at 
more usual values of y. 

With two real roots, the branch of the locus going to X ■ -Kp will go 
into the right half plane - become unstable - if Kp < 0. This root is on the 
real axis, so passes through the origin (X = 0) to go into the right half 
plane. That is, a static instability, or divergence, occurring due to the net 
spring rate being negative. The stability boundary is crossed at 
y/16 « -v 2 /Kp, and the flap motion is divergence unstable for larger y. This 
may better be viewed as a limit of Kp*. 


v 2 

Kp > - ~ (6) 

S 

required for stability. At this boundary, the other root Js at X ■ v 2 /Kp 
(hence stable, since Kp < 0) . 

8 



Hover root locations which will be of particular interest are those where 
the frequency is near & Multiple of Vrev . The hover root has a frequency of 

Vrev for y/16 ■ Kp ♦ /v 2 + K p 2 - (1/4). For Kp * 0 or small, this means 
y approximately 14, so the hover root frequency is above %/re v usually. The 
hover root has a frequency 1/rev for y/16 * Kp + /v^ ♦ K p ^ - 1. For K P = 0 
or small, this means y must be small. Since v > 1, there is only one 
crossing of ImX * 1 or V rev by the y locus (except for the case v * 1, 
when the locus starts at ImX ■ 1 for y ■ 0, and will have a second crossing 
if Kp > 0). 

The natural frequency of the flap motion is u^ 2 » v 2 + (y/8)Kp. This is 
mainly given by the structural and centrifugal stiffening, that is, by the 
frequency v; it depends on y only as it influences the effectiveness of 
pitch/flap coupling. Negative pitch/flap feedback, Kp > 0, adds a positive 
aerodynamic spring and so increases the effective flap spring rate. Notive 
that u»n * 0 gives the criterion for divergence of the flap motion. The 
damping ratio is c * (y/lb)/^. Hence, the flap motion is very heavily 
damped in hover, with c typically 50 percent critical damping, due to the 
high aerodynamic damping. 

Expansion in y 

Consider a perturbation solution for the stability in forward flight at 
small u. Using the method of multiple time scales (as described in Appen 
dix B) , the behavior of the system is examined for 4 of the order 1, y' 1 , 
y' 2 , etc.; that is, let >' 4* * “ y4» 4* ■ U 2 i>, .... Then the time 

derivative is 

JL 3 JL * JL 2 JL + 

54 34o 54; 3ip2 

Expand B as a series in y, each tern depending on all the time scales i(' n : 

6 ■ 6 0 (4o.4i* 42.- • •) + uBj (4o»^i» • ••)♦••• 
and also expand the parameter y as a series in u: 

Y " Y o * WY l * v2y 2 * * * • 

This expansion of y is a way of quantifying when y is near certain 
critical values; that is, if y 0 is some critical value, and yj ■ (y-y 0 )/v 
is order 1, then y is order y near y 0 . The quantity y is still the 
parameter given, hence this decomposition into y 0 , y 1# Y 2 » . . . changes 
with y. 

Now 6, d/d4, and y are all expanded as series in y. ’ r hese expansions 
are substituted into equation (3), and all terms of the same order in y 
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collected, assuming that all the coefficients in the expansion are of the same 
order. A fundamental assumption of the method of multiple time scales is that 
the 6 n must all be the same order for all the time scales if*.,. All the 
terms of like order in y are collected and separately set to zero, to obtain 
the equation that starts the analysis at each order. 


Order 1 Results 


The order 1 ter.r^ of equation (3) give 


» 2e o 

3 *0 2 


+ 


*0 

8 


3*0 


+ 


( y2 + K P if) 6 0 = C 


The solution of this equation is 


B 0 = Re Boi (*j , 


• ) e 




o v o 


where the root X 0 is 


\q = 



(7) 


( 8 ) 


(9) 


The convention will be followed that Xq is the root with positive frequency; 
the other root is the conjugate Xj. Then to order 1 the equation of motion, 
and so the roots, are just the hover limit. Since 8 0 is a function of all 
the time scales, equation (7) is a partial differential equation, which 
determines B 0 as a function of i|i 0 only. Thus, 6 o 1 still depends on ipj, 
<P 2 , etc. 


Order u Results 


The order y equation is 


[A 

K 2 


Y 0 3Bj 


+ T + V 2 * k p T 


) B >] 


= 2 




Y o 36 o Pi Y o . . \ 

a * 0 \T* T u sln V 

(x r js *o * X K P * k p sin *o)*o 


30r 


( 10 ) 


This is a differential equation for forced by the order 1 solution 

Bq C 4»o ) 1 substituting for B 0> it becomes 



(VBi >0 36 1 ( 2 

dlpQ 2 8 ^0 \ 



+ 



36 01 

3*1 


+ 


T 


Uc 


+ 


Kp'iRr.. 


”]■ 




0^0 


+ JJ - i(*0 * 2 Kd)J 6oi e ^° i ^° 


+ 12 f 1 + i(X ° + 2K P )] e 0 ie (Ao-l) ^° 


+ conjugate 


( 11 ) 


The right-hand side has a term of the form A]e due to the Bo solution, 

where Aj is independent of iJjq. But the left-hand side is the same as the 
order 1 equation (eq. (7)), so it has the same homogeneous solution e A °^°. 
Then the solution for 3i is of the form 


Bi 


Re 



Al 


2X 0 




+ 


( 12 ) 


The forcing of equation (10) by its own homogeneous solution produces the 
second term in 6^, which is order compared to the solution for 3 q. Su 

as <l>o increases, Bi will become arbitrarily large compared to 8 q, which 
violates the assumption that all the terms S n in the expansion of B are 
the same order. This situation can only be avoided if the secular term Ai, 
the coefficient of the homogeneous solution forcing equation (10), is 
identically zero. 

Assume for now that Xq ± i f X 0 ; then the periodic coefficients do not 
contribute to the secular term of equation (10) . Setting the secular term to 
zero gives a differential equation for Bqi(^i) : 

( 2X ° + ~i) Hr 1 + -g- oo + k p^oi = 0 cm 

or 

- moi ■ o (i4) 


where 
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(X 0 


Kp) 


A i = - 


2iImAo 


A, <(/, 

The solution of this equation is So l = 802 (^ 2 *- • .)© 1 » s ° the solution 

for 8 so far is 


g = Re jg 02 (<l>2 > • • .)e X °^ 0+Xl ^ 1 J + 0(y) 


(15) 


The eigenvalue to order y is then 


X 


Xq + uXj 



(16) 


That is, Xj is just an order y perturbation of the hover root, due to the 
order y expansion of y. To order y then, the eigenvalue remains just the 
hover value, with no influence of forward flight at all. 

The assumption that Tq ± i i Xq means that IitiXq + Vnev. This 
requirement then is that the hover root X should not have a frequency near 
Vrev; "near" means being able to write the hover root as Xq plus an incre- 
ment, such that when Xq has a frequency of Vrev, X - X 0 is order y 
small. As y increases then, the distance from Vrev which considered "near" 
increases. The analysis will return to the case ImXo = Vnev later. 


With the secular term removed, equation (10) becomes 




Bl Y 0 + (yj2 + V M n 

7 + T^ + v + k pT/ 6 i 


[l - i(X 0 + 2K p )]g 01 e (Xo+i) ' l 'O 
[l + i(X 0 ♦ 2K P )J 6 0 ie (X °~ 1 ^ 0 


+ conjugate 


(17) 


The solution of this equation is 

Bi = Re je 1 1 ('l' 1 )e X °^° ♦ Soi [A + e (X ° +i) * 0 ♦ A. e (X °' l) ’ , ’ 0 J J (18) 




\ t r 

H > ; 
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where 


1 + i (X Q + 2Kp) 
±2ImXo + 1 


V 

r * 


l 


2^! 2Tti : n! 

the solution so far is Bo^O^l) ancl 


-[• 


The order u 

86 


Order u 
2 equation is 


2 Results 


£Si. Is. 

< 8 


S' * (* 2 * 62 ] 

+ K p T + Kp T" sin lD °) 61 

sin ♦») 

. cos * 0 * ^ si" 2*0 * *P ^ sin - T C0S 2 *“)1 6 ° 


8<i»0 9 ^1 

cos ipo 

a 2 6 


a.* 2 _i!Sa_ 

2 Z 3^0 ^ 2 


(19) 


j and fii are substituted into this equation, 

2S 3T.SS iS-STf i‘ ^n^Ution. .M. 

asst«ed that thS hover frequency is not near 

4/ rev. Then the secular term is 
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where X 2 


is the order y 2 term in the expansion of the hover root 

Xq + yXi + u 2 X2 * • + i.A 2 + T Kp (lk) + °^ u ^ 


( 21 ) 
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Regarding this as a differential equation for 0 1 1 l ) > it is forced by its 
own homogeneous solution eM^l. fe ting the secular term of equation ( 20 ) 
(i.e., the entire right-hand side) w zero gives an equation for 602 (^ 2 ) '• 
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Thus, the eigenvalue to order u 2 is 
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(23) 


(24) 


*8 


Then to order y 2 - and for y and v such that the hover frequency is away 
from h/rev and 1 /rev - the influence of forward flight on the eigenvalues is 
simply a small (order y 2 ) change in the frequency. The first effect in the 
frequency just corrects the (y/ 8 )Kp spring term to account for the increase 
of the mean of KpMg with y. The second effect of y is entirely due to 
tne periodic coefficients. This order y 2 change in the frequency is in fact 
quite small for y out to 0.5 or so, as the examples below will show. There 
is no influence of forward flight at all on the damping of the root. 


This expression for the eigenvalue was derived assuming that the hover 
root, are a complex conjugate pair; it may be shown, however, that the same 
expression is valid for y large enough that the two roots are real, that is, 
when the radicand is negative. A point of particular interest is where one 
branch of the locus on the real axis crosses into the right half plane, that 
is, the divergence stability boundary. The criterion for this boundary is 
that X — 0 , which from equation (24) gives 


v 2 + (1 + y 2 ) X K p 



8 

¥ 


. I 
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Kp 


4K p 2 
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or since -(y/8)Kp is an order u 2 distance from v 2 , the criterion on Kp 
for stability is 


(1 + 



(25) 


The y effect on the right-hand side (the periodic coefficients) dominates 
that on the left-hand side (the average of KpMg) . So the critical value of 
negative K p , beyond which the locus lies in the right half plane, is actually 
increased by y; the hover criterion on Kp is then conservative in forward 
flight. That is the opposite conclusion as would have been reached consider- 
ing just the averaged coefficients. In any case, however, the influence of 
forward flight is only order y 2 small. 


Equation (24) also gives the effects of y on the boundary between where 
there are two real roots, and where the roots are cor jugate pairs. This 
boundary is given by ImX = 0, which to order y 2 is 



/v 2 + Kp 2 ) 


16 


1 + y‘ 


v 2 (-Kp + /v 2 + Kp 2 ) Kj 


/v 2 T Kn 2 


) 


(26) 


For Kp = 0 this reduces to simply 



V (* ♦ y2 ^ V 2 ) 


(27) 


Near Vrev Frequency 

Now return to_the case when the hover root frequency is near Vrev. If 
ImXo = 1/2, then Xq + i * Xq and the periodic coefficients contribute to the 
order y secular term. The criterion ImXo * 1/2 means 

T5 * ♦ /** * K P 2 - 7 • 

?ince Y = Yq + WYj + • • • » that means y must be such that Y ■ Y 0 is 
order u small. The order 1 root is then Xq * -(Yo/^) + (i/2), and the 
order y secular term is now 

^2X 0 + + (Xq + Kp)B 0 i + yy [l - i (*0 + 2K p )] 8oi a 0 
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or 


^01 

3+1 


- XjSoi 



6 P 24 / 


&01 = 0 


(28) 


where here, since ImXo = 1/2, the order u expansion of the hover root is 



The ITqj term is the periodic coefficient contribution to the secular term. 
The solution of an equation of this form is given in Appendix C. The solution 
for B q i depends on the quantity 


D 2 = U»»l ) 2 - - X Up) * (jj) ] 

Now if D 2 > 0, the solution for 8 qi has ter ® s with time behavior like 

-(ReAj±iD)+i -q (yi/16)+±iuD+ 
e * e 


and then 8 0 has terms like 


6 Vo - e " (y/16)+±i+ [( 1 / 2 )+udJ 


The damping is unchanged, and there is an order y change in the frequency, 
due to the periodic coefficients when the hover frequency is near Vrev. If 
D 2 < 0, then Bo i has terms like 

e -(ReX 1 ±D )+ 1 _ e [“P(Y 1 /16)±yD] + 
and then 80 has terms like 

BpjeVo _ e [-(Y/ 1 &)±uD]++(i/2) + 

There is an order y change in the damping, both more and less stable, while 
the frequency remains fixed at *j/rev. 
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The influence of u on the roots near Vrev frequency then is first an 
order y change in the frequency towards Vrev, with the damping at the hover 
value. When the roots reach Vrev, the frequency remains fixed while there is 
an order y change in the damping. One root is stabilized, the locus moving 
to the left on the X plane, but the other is destabilized, it moves to the 
right. This type of behavior of the roots is characteristic of periodic sys- 
tems (as discussed in Appendix A). Indeed it appears here due to the contri- 
butions of the periodic coefficients to the order y secular equation when 
the hover frequency is near Vrev. There is in this case a critical region 3 
inside of which a change in the flap damping occurs. For many problems with 
periodic coefficients, the system is unstable inside such a region. In this 
case, however, the hover damping, ReX = -y/16, is quite large; and the change 
in the damping is only ±yD, hence, order y small compared to the hover 
damping. So the critical region is a region of stability degradation rather 
than instability (for small y, i.e.). 

The boundary of the critical region is given by D 2 = 0, or 


(if ' K <>) T ■ 1 T? f 2 - T K P * V C30) 

Since y = yo + uyi> write y = yo + Ay where y is the value such that the 

hover root is at Vrev. Then if the hover frequency with y is near V rev, 

= Y - Yo roust be order y small. The boundary of the critical region 

is then 

(if " Kp ) = ±u 2T A 2 ' T k p + 4Kp2 (31) 


Considering the y root locus, then the critical region boundary is the value 
of y for which the locus reaches Vrev frequency, and is just about to 
encounter the stability change at Vrev frequency. From the behavior of the 
y root locus, this boundary value will be denoted F corner , where then 


w corner 






li 

8 


+ 4K p 2 


(32) 


For Kp = 0 this reduces to u corner ■ ±(Ay/16) (3/2v) . In general, this 

equation gives the boundary of the critical region on the y - y plane. The 
eigenvalues are 
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A = Ag + uReAj - iuD 



For y = 0, this reduces to 




' 1 fe 2 (ii * "p) ■ 


which is just an order y (i.e., order Ay) expansion of the hover root from 
Vrev frequency at Yg* 

Near 1/rev Frequency 

Consider now the case with the hover root frequency near 1/rev. With 
ImAg * 1^ then Xg + 2i = Xg, so the periodic coefficients contribute to the 
order u‘ sec ular term. T he criterion ImXg = 1 means 

y 0 /16 = Kp ♦ /v 2 Kp 2 - 1, and the order 1 root is then Xg*-(Yg/16) + i. 

For this case ImXg i 1/2, so the order y results are applicable; then the 
root to order u is known (i.e., X =* Xo + yXi), and now the orde-** y 2 influ- 
ence is sought. The secular term of equation (19) becomes, for ImXo = 1: 


(nr ' x,Su ) ■ (nr * [' X2 * (t ) 2 nr ( v2 - 

-*r T? 

* {rr (if) [ 2 • * 2Kp)] a 

- (1 - i^ p )J Bo2 eXl '* >1 


X K P + 4K P 
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The 802 term is the contribution of the periodic coefficients to the 
secular term when I 111 X 9 = 1. The secular term of this equation is the coeffi 

cient of e^ 1 ^ 1 , hence, it is exactly .$ before (eq. (24)) - giving the same 
result for the eigenvalues to order y 2 - unless \\ = Xj. In that case, 
which requires that Xj be real, the entire right-hand side is the secular 
term of equation (35), including the 802 term. Since 




requiring ImXj = 0 means that = 0; and so Xj = 0 in fact. For the 
periodic coefficients to contribute to the order y 2 secular term requires 
then that y = Yo + u 2 y 2 + • • • » that is, y must be such that Y - Yo is 
order u 2 small; this quantifies what "near" 1/rev frequency means. For a 
given u then the hover frequency must be closer to the critical value than 
was required for the Vrev region. The effect of forward flight near 1 /rev 
frequency is smaller than near Vrev, only order y 2 compared to order y 
for the latter. 
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With Yj =0, the secular term of equation (35) becomes 

* [-i re ('0 * v * (x) 2 re G 2 - ^ k p ♦ v) - k p re ‘l 602 

2 

+ [• if* iK P if ' 1 i {t) (" if* 2 k p) (' if + 2 K p + *)] 6 ° 2 


= 0 (36) 


The solution for 802 depends on the quantity D (see Appendix C) given by 


D 2 
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Gt - ^ * ( t ) 


2 v 2 - -j- K p * 4K P 2 
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p H> 


(re)([‘ * -f (re 2 k p )] 2 * [ k p • ^ (If - ihj ] } on 


The behavior of the roots is similar to that near Vrev, except that all 
changes due to forward flight are here order y 2 . There is a critical region, 
inside which there is an order y 2 change in the damping of the flap motion. 
For small u, the damping is the same as the hover root while there is an 
order y 2 change in the frequency, towards 1/rev. At the boundary of the 
critical region the root reaches 1 /rev frequency; and for still larger y, 
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inside the critical region, the frequency is fixed at 1/rev while there is an 
order j 2 change in the damping, one root becoming more stable and the other 
less. 

The boundary of the critical region is given by D 2 = 0. Since 
y * y 0 ♦ u 2 y 2 + . . ., write y = yq + Ay where y 0 is the value such that 
the hover root is at 1/rev. Then, if the hover frequency with y is near 
1/rev, iv = i - yq must be order u 2 small. The critical region boundary 
may then be written 

(fe ' K p) 16 = p2(Cl 1 C2) (38) 


where 


C, * 


&r- 


- -f Kp ♦ 4 K f 
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+ Ki 
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3? «/[* + T (l& ■ 2 k p)] * [ K P ' TT GT 2K p) ] 

J\ l + T ^ ■ l) ] + [ Kp " T (v2 ' 1} (li " 2 K p)] (40) 


Then the u root locus reaches the critical region boundary at 
^corner = - ! or ^ where 




(is ■ k p) t£ 


Ci ± C 2 


(41) 


For K-i * 0, this becomes 
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ipj v 2 /v 2 - 1 ± J /'l'* (v 2 - 1) ♦ v 2 (\> 2 - l) 2 


(42) 
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The eigenvalues are now 


X = x o "*■ y 2 ReX 2 - iy 2 D 
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k * i * ‘* 2 


16 


+ 1 


/[rf (rk - k p) - c >] 2 - c * 2 

1 16 (l6 ' K p) /(/ ‘ C 7 ) (‘ ' 


(43) 


For u = 0, this reduces to simply an order Ai|» (order u 2 ) expansion of the 
hover root f’-om 1/rev at y & . 

Summary, and Discussion of the Results 

This section summarizes the results of the perturbation solution for the 
influence of forward flight on the helicopter rotor blade flapping stability. 
The expressions for the eigenvalues are collected from the analysis above. 

The behavior of the roots is discussed, in terms of the y root loci and the 
critical regions on the y - y plane. These results are for the shaft-fixed 
stability of an individual rotor blade. 


Hover.- The hover limit, y = 0, has the eigenvalues 



(44) 




i 


These roots are usually a comp lex con jugate pair, located at ReX = -y/16 on 

the circular arc with radius /v 2 ♦ Kp 2 and center at X = -Kp. Forward 
flight, u » 0, introduces periodic aerodynamic forces into the dynamics, which 
radically influences the behavior of the eigenvalues and the analysis tech- 
nique required to obtain them. A perturbation method based on small y has 
been used to obtain explicit expressions for the roots when y > 0, including 
the effects of the periodic coefficients. It is an order y 2 analysis 
(consistent with the neglect of the reverse flow region effects), which is 
valid to approximately y ■ 0.5. 


Forward flight, away from critical regions.- When the hover root 
frequency is not too close to a multiple of Vrev, the roots to order y 2 ire 


X 
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v 2 - J Kp ♦ 4K P 2 
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(45) 
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There is only an order u 2 change in the frequency, which is quite small even 
up to u = 0.5. This expression applies in particular when there are two real 
roots, that is, when y is large enough so the radicand is negative, so it 
gives the criterion for divergence instability, including the influence of 
forward flight. The divergence boundary is given by where one branch of the 
locus on the real axis goes through the origin into the right half plane. The 
boundar> criterion is then X = 0, for which equation (45) gives the criterion 
for divergence stability as 


(1 ♦ U 2 ) ^ K P > 


-v 


1 ♦ c 



(46) 


This is a limit on the pitch/flap feedback allowed; the divergence instability 
occurs when Kp is sufficiently large negative (positive feedback). A con- 
stant coefficient approximation to the equation of motion, using the average 
of the coefficients, includes only the effect of the mean of M 0 on the left- 
hand side of equation (46); that is, a conservative approximation for this 
case. 


Forward flight , near h/rev frequency.- It is characteristic of a system 
with periodic coefficients that for certain values of the parameters there 
occurs a degradation of the stability. Typically this occurs where the basic 
eigenvalue - here the hover root - has a frequency corresponding to a multiple 
of one half the fundamental frequency of the equation coefficients. Specifi- 
cally, consider when the frequency of the hover root of the flap motion is 
near Vrev. There occurs then an order u influence of forward flight. When 
the frequency is near V*ev, equation (45) is no longer valid r or X, and the 
following result must be used instead. Let y 0 be the value of y for which 
the hover root would (with the given v and Kp) have a frequency exactly Vrev; 

that is, Yq/16 » Kp + A > 2 + Kp 2 - (1/4). Write Ay * y - y^. Then, if 
Ay/ 16 is ordei u small, the eigenvalue is given by 


& * r - 1 ft 2 ( t £ - "p )/ 1 - 


(47) 


where 


& (if - ^ 

y comer * Y 

if A 2 - T Kp * 4lt p J 


(48) 


The subscript "comer" refers to the behavior of the v root locus on the X 
plane, as discussed below; it is the boun dary o f the critical region. For 
Kp ■ 0 this result reduces to v 0 /16 = - 1/4 , and u corn er * (Ay/16) (3/2v) 

for the boundary. 
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This solution exhibits the .v*ii:g behavior. If u < ,-ner» there is 
■> order u change in the freer:-. - , toward V rev, while the real part of the 
root remains fixed at the » iae of -y/16. At y 3 Ucorner the roots 

reach r m • - Vrev. For *ner the frea ency remains fixed at 4/rev 

while the.t- is an order >. - .s’»ge in the damping; ReX is increased for one 

root and decreased for the c ». ( ev. This behavior of the locus near 
ImX - Vrev is character 1st if, of periodic systems (Appendix A). While there 
is a stability degradation if u is large enough (greater than Ucorner* 
which decreases with A'., i.c., as the hover root approaches Vrev), the 
reduction in damping is order u small. The hover damping, ReX = -y/16, is 
quite large for usual values of y* and so stability is maintained for small 
u, even with the influence of the periodic coefficients. 


Forward flight , near 1/rev frequency.- Similar behavior is exhibited when 
the hover root frequency is near 1/rev. Let Ay = y - Yo» where now y 0 is 
the value of y for which the hover root (with the g iven v and Kp) is 
exactly at 1/rev, that is, y 0 /16 3 Kp + /v z + Kp 2 - 1. Then if Ay/16 is 
order small the roots are given by 



where the corner u are 


2 2 

U1 ,U2 = 



Cl 1 c 2 


(SO) 
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expressions for Ci and C 2 arc given above (eqs. (39) and (40)), along with 
the results for the limit Kp 3 0. The behavior of the u loci is like that 
near Vrev, except that here all changes are only order u 2 . For small u 
there is an order u 2 change in the frequency, toward 1/rev, while there is 
no cl ange in the real part from the hover value. For u ■ ui or u 2 (only one 
will be real) the roots rea^h ImX = 1/rev, the critical region boundary. For 
still larger u, the frequency is fixed at 1/rev while there is an increase of 
the real part of one root and a decrease of the other. The stability degrada- 
tion is only order u 2 small, so again the flap motion will remain stable for 
small u and reasonable y. 

Comparisons with numerical solutions for the flap roots (from calcula- 
tions by the author, and from the literature, e.g., ref. 12) indicate that the 
perturbation solution to order u 2 is accurate to about u * 0.5. This 
solution then covers the range of interest for most helicopters. 

Root loci for forward flight.- Typical y root loci are shown in 
figure 1, for several cases of v and y. The cases considered are: (a) v ■ 1 
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b 


i , 


and y « 10; (b) v * 1.1 and y * 6; and (c) v ■ 1 and y » 6. The pitch/flap 

feedback Kp => 0 for all three cases. Then the hover roofs, y 0, are 

located on a circle with radius v and center at the origin. Case (a) is a 
typical articulated blade; with the large y and v ■ 1, the hover frequency is 

well below 1/rev. Hence, with y > 0 the locus encounters the 4/rev critical 

region. Case (b) is a typical cantilever blade (i.e., a hingeless rotor); 
with v » 1 and a lower y, the hover frequency is above 1/re/. Hence, the 
locus encounters the 1/rev critical region. Case (c) is an example « s f the 
behavior when the hover frequency is away from any multiple of ^/rev. So the 
only influence of y is a very small - order y 2 - change in the frequency. 
The loci shown in figure 1 thus cover all the cases of away from the critical 
regions, near tyrev, and near 1/rev. They illustrate the form and magnitude 
of the influence of forward flight. Specifically, th^ behavior in the 
critical regions at V rev and 1/rev frequencies is shorn; and the large - 
order u - effect cf the 4/ rev region, and the quite small effect away from all 
critical regions may be seen. 

y - y plane.- The eigenvalues depend primarily on y and y, so the above 
results may be presented as contours of constant ReX and the constant ImX 
on the y - y plane. This is the presentation found in much of the litera- 
ture. Such plots are shown in figures 2 to 5 for Kp = 0 and v ■ 1, 1.05, 

1.1, and 1.15, respectively; and for v * 1 and Kp = 0.1 in figure 6. They 
are based on the perturbation solution given here. The expressions for the 
critical region boundaries in terms of the comer y for thj y root loci 
may be rearranged, to give the boundary in terms of y for a given value of 
u. Writing Y 3 Yq + Ay, the two real root region boundary is 



where y 0 /16 » Kp ♦ /v 2 ->■ Kp 2 . The *s/rev cr tical region boundary is 



(51) 


(S2) 


where y q /16 * K p ♦ /v 2 ♦ K p 2 - 1/4. Ani the )/rev critical region boundary 
is 



t 
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(53) 


where > 0 /16 ■ K p ♦ /v 2 ♦ Kp 2 - 1, and the constants Ci aim C 2 are given 
above (eqs. (39) and (40)). This expression shows that Ci gives the offset 
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of the 1/rev critical region boundary frost y and C 2 gives the width of 
the region; the Vrev region has no offset to°ordor y, so is symmetrical in 
Ay about y Q . For Kp = 0 these boundaries reduce to 



with = /v 2 - 1 (56) 


respectively. The boundaries in figures 2 to 6 were constructed using these 
expressions Only the critical region and real root boundaries are shown 
here, but expressions for constant ReA and ImA may also be obtained from the 
perturbation solutions. The presentations of the y - y plane results given 
in the literature usually include the ReA contours at least (since only 
numerical results are usually available, the y - y plane is the most 
efficient way to indicate the stability trends with y and y variations). 

No instabilities are encountered on the y - y plane, for the range of 
parameters shown. A divergence (static) instability is encountered at hign 
y if Kp < 0, and an instability in the critical region (usually the 1/rev 
region first) is encountered for much higher y (around 2). The critical 
regions, combinations of y and u where the frequency is fixed at \ or 1/rev, 
are the principle effect of the periodic coefficients. Notice they encompass 
more and more of the y range as y increases, that is, as the periodic 
coefficients increase. The two real root region (ImA = 0, i.e., two roots on 
the ReA axis) is due to y being large enough so the flap motion has super- 
critical damping; while it is influenced by the periodic coefficients, it is 
not the same type of phenomenon as the critical regions (the roots for this 
case are given by the same expression as for two complex roots away from the 
critical regions). 

A horizontal line on the y - y plane is a line of constant y, hence, 
the variation of ReA and ImA as such a line is traversed gives the coot 
locus for varying y. For example, consider the y locus for case (a) above, 
= 1 and y = 10; the y = 10 line is indicated on figure 2. As y increases 
from zero, the line remains parallel to ReA = constant lines, so ReA 
remains fixed at the hover value. The ImA * Vrev critical region comes 
closer to the horizontal line, indicating that the frequency of the root 
approaches Vrev. At the corner y, the locus crosses into the Vrev region. 
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For higher y the locus is in the critical region, so the frequency reaains 
fixed at Vrev while for each point in the region there will be two values of 
ReA, one more and one less stable than the hover root. The constant y lines 
for cases (b) and (c) are also shown in figures 2 and 4; the behavior of the 
y - u plane along these lines ma> be compared to the u loci shown in 
figure 1. 

Figure 7 shows the y - y pUne for v = 1 and Kp = 0, from a 
perturbation solution equivalent to the present c-ne but including the order 
u 2 correction to the Vrev critical region boundary. That solution, from 
reference b, gives the Vrev boundary as 


Ay 
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(55a) 


to order u 2 . Figure 7 may be compared to figure 2, which is based on the 
present solution, hence, only to order y in the Vrev region. While the 
order u 2 influence is not negligible for y = 0.5 or so, the major effects 
of forward flight are contained in the order y solution. The results of 
reference 6 are discussed further below. 

High y behavior.- For y larger than 0.5 or so, a numerical method 
must be used to calculate the eigenvalues of the flap motion. It is also 
necessary then to include the reverse flow effect in the aerodynamic flap 
moments. At very high y (above 3 to 4, say) a perturbation solution based on 
an expansion in y' 1 is possible (ref. 14). Such a solution is of less use 
than the small y solution, at least for current helicopters; it does, how- 
ever, give some insight into the high y behavior (refs. 14 and 15). For y 
around 1, only numerical solutions are possible, unless some other parameter 
is used for the perturbation variable (such as y, as in ref. 14). 

At u = 2.2 or so - the exact value depends on v, y, and Kp, but there 
is not much variation for the range of parameters of current helicopters - a 
flapping instability is encountered. It occurs in the 1/rev critical region, 
the root being destabilized by the periodic coefficient influence crossing the 
imaginary axis into the right half plane. The y - y plane for v = 1 and 
Kp = 0, with u out to 2.5, is shown in figure 8; this plot is a composite of 
the results available in the literature. Figures 2 and 7 give the y - y 
plane for the same v and Kp, but for y to 0.5 only. The behavior of the 
critical region boundaries at high y, and the high y instability in the 
1/rev region are shown. The high y behavior of the y root loci in 
figure 1 may also be inferred from figure 8. 

For case (c) , v ■ 1.1 and y * 6, the y locus enters the 1/rev critical 
region at small y. The locus remains in that region, one root moving to the 
left and the other to the right, until the latter branch becomes unstable 
(crosses the imaginary axis into the right half plane) at high y. Fo* 
case ib) , * * 1 and y » 6, the locus shows a decrease in frequency up to 
u ■ 0.5. But at somewhat higher y the loci turn around and the frequency 
begins to increase, for at about y ■ 1.1 the roots encounter the 1/rev 
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critical region - see the y = 6 line in figure 8. In the 1/rev critical 
region for higher p , one branch of the locus is destabilized, until an 
instability is encountered at about p = 2.35. 

For case (a), v = 1 and y = 10, the u loci also encounter the high p 
instability in the 1/rev critical region, so first the roots must get from the 
Vrev region to the 1/rev region (see the y = 10 line in figure 8). As p 
increases above 0.5, the loci in the Vrev region eventually turn around, the 
real parts of the two branches then approaching each other instead of diverg- 
ing. At p about 1.55 the roots get back to the hover value of damping, and 
then break away from the Vrev region. The frequency of the roots increases 
toward 1/rev then, while the damping is fixed; that is, the roots are complex 
conjugates again. When the roots reach ./rev they enter the 1/rev critical 
region. The transition from Vrev to 1/rev frequency occurs very quickly, 
during a very small p increase, because the corridor between the two criti- 
cal regions is very narrow at this point (fig. 8). In the 1/rev critical 
region for higher p, one branch of the locus is destabilized then, and 
eventually a flapping instability is encountered at about u * 2.25. 


N-BLADED ROTOR EQUATIONS OF MOTION 


Consider a rotor with N independent blades, with no coupling by rbaft 
motion and only the excitation due to the blade pitch control. The flap 
motion of the rotor is described by a set of N equations, each of the form 
of equation (3): 


g« + (| + £ u sin ^)e (m) + jv 2 + p ccs + } u sin *m) 


+ K P [s + 3 u sin * 4 (V sin ’M 2 ]) 00 ” 3 


* [ 8 + i u sin + 4 sin C 57 ) 

where 8^ is the flap degree of freedom for the mth blade, m ■ 1, . . ,,,N. 
The azimuth location of the mth blade is \J> m « tp + mAij>, AiJ) = 2ir/N. 

For a teetering or gimballed rotor, the blades do not act independently, 
so the rotor motion is not described by these equations. They may be used, 
however, to derive the appropriate equations of motion. 
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TEETERING ROTOR 


Equation of Motion 

Consider a teetering rotor: a two-bladed rotor with a single flap hinge 
at the center of rotation. The two blades are not independent then, rather 
the rotor flaps as a whole, one blade up and one down. The equation of motion 
must be obtained from equilibrium of moments on the entire rotor rather than 
on individual blades. The coning motion of the rotor - both blades up or both 
down at the same time - is reacted by the structural restraint at the blade 
root. That is, it is a cantilever mode type of motion, for a very stiff 
blade, so with a very high frequency. In the teetering motion, however, the 
blade acts like a hinged rotor, usually, in fact, with no hub spring at all so 
the flap natural frequency v = 1/rev (although the general case of v >. 1 
will be considered). Consequently, the coning motion of the rotor will be 
neglected as a higher frequency motion, and only the teetering degree of 
freedom considered. 

Let B be the degree of freedom for the rotor teetering motion, so 
8^ = B and 6^ = -B. Equilibrium of flap moments on the entire rotor is 

given by half the difference between the B*- 2 ) and B^ 1 ) equations of motion. 
Then, the equation of motion for the teetering rotor flap motion is 


8 ♦ J 8 + [v 2 + J n 2 sin 2 <|> + Kp(I + £ y 2 - J y 2 cos 2 ^)]b - 0 (58) 


Notice that all 1/rev harmonics drop from the coefficients of 8, because 
with the teetering motion the rotor has a period of only it with respect to 
the aerodynamic environment. As a consequence, all the order y terms have 
dropped, leaving only the order y 2 influence of forward flight on the 
coefficients of the flap equation. 


Hover 


For the hover limit y = 0, the roots are as before (for independent 
blades) : 


X = .. JL 

16 




(59) 
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Only order 
is used now: 


This is the 


with solution 


where the order 


The order 


Expansion in y 2 

y 2 terms appear in the flap equation, so an expansion in y 2 


_3_ 

3* 



S = 8q + U 2 &2 ♦ • • • 


Y = Y 0 ♦ U z Y 2 + • • 

Order 1 Results 
hover limit again: 


r!o 

^o 2 


^ ^ + ( v2 + Kp **■) B ° = 0 


$0 * Re[e 0 2C'l'2»- • 


L eigenvalue is 


16 


o . . / ..2 


+ 1 / v‘ + 


la . M 

T K f ' \rsj 


Order y 2 Results 
j 2 equation is 


(60) 


(61) 


(62) 
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Yq Wz t 
8 3i|>q 



^ Yo 9&o + Yj_ 98 0 
3^0 ^2 8 3^2 8 3i/)q 


C* 


sin 2*o ♦ K p ^ * Kp ~ - Kp cos 2* 0 ) 0 O 




* C-i - KpjSc^O* 21 )*., 

♦ ^5 (i - ^p)6o2®^^~^ 1 ^ 0 


+ conjugate (63) 

Assuming that Xq + 2i t Xq, that is, IuXq i 1/rev, then the secular term is 

(2X 0 * ^f) ♦ [x Uo + Kp) + ^ K p ]b 0 2 s 0 

or 


30 

- x 2 e 0 2 ■ 0 C 64 ) 

where 

-j- (Xo + Kp) + -y- Kp 
X? - * ' 2iImX 0 

The solution is 

602 » 604C«4.- • .)e Xa * 2 (65) 


J 


V 

J 


J 


* 



* - * • 


* 


£■ 


30 



So 


6 = Re(e 04 e X °^ +x 2 'l' 2 ) + 0 (y 2 ) 


( 66 ) 


Then, for ImXo i 1, the eigenvalue to order y 2 is 

X = Xq + y 2 X2 




(1 + P 2 ) l K p 



(67) 


The only influence of forward flight on the eigenvalues of the teetering rotor 
is then the Increase in the mean of KpMg (to order y 2 , and for ImXg 1). 
This result is substantially simpler than for the independent blade, with much 
less influence of the periodic coefficients. The difference is the result of 
the internal cancelling of the 1/rev, order y flap moments for the teetering 
rotor. In fact, if Kp = 0 there is no order y 2 influence on the roots at 
all, they remain at the hover values; if Kp > 0 there is an order y 2 
increase in the frequency, due to the increased effectiveness of Kp acting 
through Mq . 

The divergence criterion (X = 0) , becomes now, for stability 

(1 ♦ U 2 ) £ K p > -v 2 (68) 


This is a stricter requirement on Kp than for the independent blade (i.e., 
the boundary is reached at a negative Kp of smaller magnitude), but the 
difference is only order y 2 . The two real root boundary (ImX = 0) is 


JL = 

16 




Kp \ 


(69) 


or for Kp = 0, y/16 a v. 


Near 1/rev Frequency 

Return now to the case Tq + 2i ■ X 0 , that i s, ImXp » 1, the hover 

frequency near 1/rev. This means yo/16 * Kp + /v 2 ♦ Kp 2 - 1, and the order 1 
root is X p * -(yq/ 16) + i. Now the periodic coefficients contribute to the 
secular term, to give 
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IjT ‘ A26 ° 2 + 1 3l (i + K P )002 = ° 


(70) 


where since ImAo 3 1, 


X 2 3 




The solution depends on 


D 2 





(71) 


The critical region boundar) is given by D 2 = 0, or writing y 3 Yg + 
with Ay order y 2 small: 


/rr ^) 


(72) 


In terms of the y locus, this gives the comer y: 


2 2 
VI >U2 




(73) 


- 3 * 3 ^ 


This expression is similar to that for the 1/rev region boundary of the 
independent blade (eq. (41)), but with 


Cl » *p if I 
c 2 ■ 37 ^ ♦ K P 2 j 


(74) 


that is, with considerably less influence of the periodic coefficients. The 
eigenvalue is 


X = Xq + y 2 ReX 2 - iy 2 D 



+ i 


Ay 

16 


' Kp )/( l ' "' 2 )( 1 ‘ ” 2 2 ) 


(75) 


which is identical to the independent blade result (eq. (43)), except that the 
comer u have diffe rent va lues. For Kp = 0, the criterion ImXp = 1 

reduces to y q /16 - X> 2 - 1, the 1/rev boundary is yi 2 ,U2 Z = ±2(Ay/16), and 
the eigenvalue reduces to 



Fot v * 1 (and Kp = 0) this is the same boundary as for the independent 
blade. For « > 1, the critical region has a smaller width, and is not offset 
(Ci = 0 still, since Kp = 0). 


Summary 

The influence of forward flight on the teetering rotor differs from that 
on the individual blade principally in that: there is no ^/rev critical 

region at all; and the periodic coefficient influence is much simpler. Away 
from the 1/rev region, the roots are given by 

x . - x . i yvToT | Kp - (k) 2 C77) 

including the forward flight influence to order u 2 * The corresponding 
criterion for divergence stability is 

(1 ♦ u 2 ) J K p > -v 2 (78) 


There is no Vrev critical region, so this expression holds also when the 
hover frequency is near Vrev. A critical region is encountered only if the 
hover frequency is near 1/rev. That is, if Ay * Y - Yo i- s order u 2 small, 

where v 0 /16 * Kp + /v 2 ♦ Kp 2 - 1, then the roots are 




♦ 1 




(79) 
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where the comer y are 


Pl 2 ,U2 2 


Ay/Io k \ 

= 16 \16 Kp y 

C 1 i C 2 


(80) 


and the constants Ci and C 2 are given above t,eq. (74)). For K p = 0, this 
reduces to u 2 omer = ± 2(Ay/16). 

The y - y plane boundaries are, writing y = y 0 + Ay, for the two real 
root region 


Ay _ 2 Y 0 P 


( 81 ) 


where y 0 /16 = Kp + /v 2 + Kp 2 ; and for the 1/rev critical region 

Ay 9 
* 




(82) 


_ 0 . 

16 


rr - K P 


where y 0 /16 = Kp /v 2 ♦ Kp 2 - 1. For Kp = 0 these boundaries reduce to 


Ay 

16 


. , Yq 
with jg * v 


(83) 


= ip 2 with * /v 2 - 1 (84) 

respectively. These results may be compart 1 with the y - y plane boundaries 
of the independent blade, as illustrated in figures 2 to 6. For Kp = 0, the 
two real root boundary is now a constant y line, that is, a horizontal line; 
there is no Vrev region at all; and the 1/rev region is the same when v ■ 1, 
but narrower and not offset (C; = 0) when v > 1. 


Consider the u root loci of the teetering rotor. When Kp * 0, as in 
figure 1 for the independent blade, then for case (a) with v » 1 and y * 10, 
and case (b) with v = 1 and y « 6 - that is, away from the 1/rev critical 
region - there is no effect of u at all to order y 2 . The locus remains at 
the hover value. For case (c) with v * 1.1 and y * 6, that is, near 1/rev, 
the behavior is basically the same as for the independent blade. However, the 
critical region is narrower now, and not offset, resulting in a higher value 
for w corner (0.407 compared to 0.286 for the independent blade). Notice that 
if the hover root frequency were below 1/rev, then Mcomer probably 

decrease because of the absence of the critical region offset with the 
teetering rotor; this is best seen on the y - y plane, as in figure 4. 
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There is then considerably less influence of the periodic aerodynamic 
forces on the teetering rotor flap stability, as compared to the independent 
blade. This is the result of the internal cancelling of all the 1/rev, order 
u flap moments. Specifically, there is no Vrev critical region, because of 
the absence of 1/rev harmonics in the coefficients; the 1/rev critical region 
is narrower, in general, and not offset if Kp = 0; and there is less influence 
of u on the two real root boundary. If Kp = 0, the two real root boundary 
does not change at all with u; this corresponds to the absence of a %/rev 
critical region, which for the independent blade pushes the 0/rev boundary up 
(as in figure 2) . 


In fact, unless the frequency of the hover root is near 1/rev - and the 
1/rev critical region is quite narrow, being only order u 2 wide - then using 
a constant coefficient approximation gives exactly the correct eigenvalues to 
order u 2 . That is, if the periodic coefficients in equation (S8) are simply 
dropped - and they are the reason for all this perturbation analysis - then 
the only effect of u retained is the increase of the mean of Mg by the 
factor (1 + u 2 ); but that is the only effect that appears in the correct roots 
for the teetering lotor anyway. This suggests that such a constant coeffi- 
cient approximation may be an adequate representation for the teetering rotor 
flapping dynamics (if the frequency of the root is kept away from 1/rev). 


GIMBALLED ROTOR, THREE BLADES 


I 



1 


Equations of Motion 

Consider a gimballed rotor: three or more blades attached to a hub with 

cantilever root restraint, and the hub to the rotor shaft by a universal 
joint. The possibility of a hub spring is included, so v > 1 is allowed 
still. The blades do not move independently now, rather the entire rotor 
moves as a whole about the gimbal bearings. There are two degrees of freedom 
describing the rotor flap motion: longitudinal tip path plane tilt Bi c and 

lateral tip path plane tilt Bis* The variables describe the rotor motion in 
the nonrotating frame; B lc is defined positive for tilt forward, and Bis is 
positive for tilt toward the retreating side. All other modes of motion of 
the rotor blades - such as the coning motion - are reacted by the cantilever 
root stiffness, so will be neglected as higher frequency motions. 

The flap motion of the mth blade is given now by 


B 


(m) 


Blc cos ♦m ♦ Bis sin ^ 


( 85 ) 


The equations of motion for the gimballed rotor tip path plane tilt are 
1 obtained from equilibrium of the pitch and roll moments on the hub. The equa- 

>. tions may be obtained from the flap moment equilibrium for the mth blade 

i (eq. (57)) by the operators 

/ 
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( 86 ) 


N 2 N 

77 2^ (• • .)cos jjj* 23 (• • *) s i n 

m=i m= l 


The result of the summation operation depends on N when, as here, the equa- 
tions have periodic coefficients. For a three bladed gimballed rotor, N = 3, 
the result is 



^ + u ^ Sin 3* 
-2 - u cos Zip 


2 - w e* ' 3i(i 


I - u X S in 3* 




The cases of N * 4 and N 2 , 5 will be considered in later sections. The 
notation 



r> 


will be used for the degrees of freedom. 


Hover 


In the hover limit, u ■ 0, the differential equation is 



i 2 ' 

• 


i ♦ 

6 ♦ 


: 2 t 


. 1 V 2 - 1 ♦ Kp J 


( 88 ) 
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The characteristic equation is 


(x 2 ♦ ^ X ♦ v 2 - 1 ♦ Kp ^ + ^2X * ^ "0 (89) 

which has the solutions X * Xq ± i and their conjugates, where 


X 0 ..^ +i /v 2 >J Kp -(^) 2 (90) 

is the rotating hover root. So the hover roots for the gimbal degrees of 
freedom, which are in the nonrotating frame, are just the rotating roots 
±l/rev frequency due to the transfer from rotating to nonrotating coordinates. 
This hover result XfjR * X R ± i is true for any number of blades, in fact. In 
forward flight the result is not this simple because the blades of a gimbal led 
rotor are not independent. 

There are two degrees of freedom, hence, a total of four eigenvalues. 
There is a high frequency mode with roots Xq + i and its conjugate, that is, 
at frequency ImXo ♦ 1/rev; and a low frequency mode with roots Xo i and 
its conjugate, that is, at frequency ImXo - 1/rev. The correspondence 
between the rotating and nonrotating regions is then as follows: 

Rotating frequency Nonrotating frequency 

0/rev 1/rev 

Vrev h and 3/ 2/rev 

1/rev 0 and 2/rev 

The 0/rev rotating frequency means two real roots; it includes the case of 
divergence instability. The regions, boundaries, and behavior of the roots 
will be discussed in terms of the rotating frequencies, so that the results 
may be compared with those of the single blade analysis. 

The eigenvectors for the roots Xq ± i are 



that is, $ic/6is B± i! this corresponds to a wobbling or whirling notion of 
the tip path plane. 


Expansion in u 


Let 
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Order 1 Results 

The order 1 equation is just the hover limit as usual: 

6 0 • 0 ( 92 ) 


The order 1 eigenvalues are then X ■ Xq ± i and the conjugates, where 



is the order 1 rotating hover root. The corresponding eigenvectors are 
S lc /8is ■ ±i, so the solution for Bo is 

J. CM, 

The first term (subscript ♦) is the high frequency mode, the second (subscript 
-) the low frequency mode; Boi± ar ® complex functions of the higher time 
scales etc - 



Ordev u Results 


The order u equation is 
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-i 


Y„ 



(Xq - i + 2Kp)Sgi+ 



e (X 0 -2i)^ 0 


+ [( 2x ° + if) ivr + ^8* (x ° + k p )s oi-] ^ J e(x °’ 1) *° 

+i yj (X o + i + 2Kp)8oi- C) e (X 0 + 2 i)i//o 


+ conjugate 


(95) 


The secular terras are the coefficients of the homogeneous solution 

^ e CX 0 ±iD^o 

on the right-hand side of equation (95) . Notice there are two homogeneous 
solutions now (really four, including the complex conjugates), as in equa- 
tion (94), and the coefficients of each must be set to zero. The resulting 
equations give the order u corrections to the roots for the high and low 
frequency modes. 
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Assuming for now that ImXo f Vrev, then the secular terms are 



or 
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01 + 




- XiB 


l p 01± 


(96) 


where 


(*o + K p) 

*1 » - — 

2xImX 0 

The solution is 8 qi± 55 8o2±(^2»‘ • Oe* 1 ^ 1 , so the solution for 8 q is 


8q = Re 




(97) 


Hence, the eigenvalues to order p are 


X = Xq ± i + pXj 


iV + 1 


p- ' \ *» • (*/ * 1 


(98) 


which is just the rotating hover root ±i still. With the secular terms set 
to zero, the solution of equation (95) is 



( 99 ) 


where 


Y 0 (X 0 ♦ 2Kp % i)(*i) 
*±““17 ±2ImX 0 - 1 

which completes the solution to order p. 
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Order y 2 Results 


The order y 2 equation is 
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Y, 
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cos 3*o + K p — sin 3* 0 sin 3 * 0 - Kp ~ cos 3* 0 


sin 3* 0 - Kp cos 3*o - cos 3* 0 - Kp ^ sin 3* 0 J 
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|(- 4 ) 
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Y, Y 0 

Y (*0 + K p) + ~g 


(». 4) 

Kp + (^0 + Z^iA+j^oi+jl^j* 


(Xo + i)*o 


. g U • K p ) Bo u] Qe (> "* l)n 

+ (-i)(X 0 + 2Kp-i)6 11+ +[^ (-i) + (2X 0 +^- 2 i)A + ]- 5 ~ i 

+ (-i)(X 0 + 2Kp - i) +-j- (Xq + Kp - i)A + Jboi+| ^ ^ ® 
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(X 0 + Kp)8u- + 



Moi- , a2g oi- 

3*2 a*! 2 


ex, * *, (.«*.] N,.) (■;) 

*[lt (-i*K p )8 0 i-]^ 1 J e(l °' i) *° 
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Assuming that ImX 0 i 1, the secular tem is 
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3\|»1 
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l p ll± 
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*) 


02 ± 


3^o 


- A o + 


f £ Kp + ^ (X 0 ♦ 2Kp) (±i)A+ 


2iImX 0 


0O2±?e AlV ' 1 (101) 


and the secular term of this is 


_ 02 ± 

3^2 


-X 2 


"8* Kp + TI + 2Kp)(±i)A+ 

2iImXg 


8q2± = 0 


It follows that the eigenvalues to order p 2 are 

X K P + tI (*0 + 2K p ) (±i)A ± 


= ±i + X 0 + uX 1 + u 2 |^ 2 + i 

- ±i-y£+i Jv 2 + (1 ♦ p 2 ) X Kp - + u 2 i 


I| (X 0 ♦ 2Kp) (±i)A ± 
2ImX n 


Which may be reduced to 


X - ±i 


- 2 (*) 2 


2 - 


v2 "8 K P + 4K P 2 


L 

v 2 + (1 + p 2 ) Kp - ^ 

X , 

±iy2 (lfr) fe ' 2K P V 


16 +i 


1/2 


1 + 


j / 

2 4 v2 - 8 K P + 4Kp2 1 


( 102 ) 


(103) 


(104) 


and the conjugates. 

Although derived here on the basis of complex roots, it may be 
demonstrated that this expression is also valid when y is large enough so 
that the rotating hover roots lie on the real axis. This single composite 
expression then is valid A r all cases except when the hover root frequency is 
near 4 or 1/rev. This result has some of the effects that have been seen for 
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the independent Made: the factor (1 + y 2 ) correcting for the increase in the 

mean of KpMe appears as usual, and the second order y 2 term in the radicand 
is like the corresponding term for the independent blade (eq. (24)), although 
half the magnitude. There are also now some order y 2 effects that give 
quite new behavior, however. There is an order y 2 change in the ±i shift of 
the nonrotating roots. Secondly, the last order y 2 term in the radicand 
produces, for the case of complex roots (i.e., as written), an order y 2 
change in the damping: an increase in the damping of the high frequency mode, 
and a decrease for the low frequency mode. This is quite different from the 
behavior seen so far, where the damping has always remained fixed at the hover 
value -y/16 when the roots are outside the critical regions. For the case of 
the rotating roots being real (i.e., bringing the i inside the radicand) 
there is an order y 2 change in the frequency of the roots: the more stable 
root going to highei frequency, and the less stable root to lower frequency. 
These new effects will be seen again when the constant coefficient approxima- 
tion is considered; a discussion ol their origin will be put off to that time, 
since that is a problem of more general interest. The present discussion will 
concentrate on outline the behavior of the eigenvalues produced by these new 
terms. 

Consider the order y 2 influence of forward flight on the divergence 
boundary. When the rotating hover roots are real - the divergence boundary 
being a special case of that - the last term in the radicand of equation (104) 
produces an order y 2 change in the frequency. Hence 




v 2 - J Kp + 4Kp 2 


1/2 


V 2 - (1 + U 2 ) J K p j 


Setting ReX * 0 for the divergence boundary, there follows the criterion for 
divergence stability: 


(1 + y 2 ) X K p > 



(105) 


The order y 2 influence on the right-hand side is just half the result for 
the independent blade, so the criterion on negative Kp is more strict in 
this case. The difference is only order y 2 , however. In the nonrotating 
frame, this divergence instability occurs at a frequency order y 2 above 
1/rev. 

Consider now the behavior of the roots of the three-bladed gimballed 
rotor near the boundary for two real roots (rotating) . Recall the result for 
the hover roots of individual blades, in the rotating frame; fig ure 9(a) 
illustrates the behavior, for variations in y. For y/16 = Kp + + Kp2 # the 

two roots meet at the real axis, and then for larger y they proceed in 
opposite directions along the real axis. The transformation to the 
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nonrotating frame simply shifts this behavior by ±i, so that the above 
behavior occurs at ±l/rev instead of on the real axis. The hover roots of 
individual blades, will in the nonrotating frame thus have the behavior shown 
in figure 9(b). This is also the behavior of the hover roots of the three- 
bladed gimballed rotor, for y near the rotating real root boundary. This 
should not be confused with the behavior of roots near a critical region; it 
occurs here at ±l/rev because of the transformation to the nonrotating frame. 
The critical region behavior involves two roots, one at positive frequency and 
one at negative frequency, which after crossing the critical region boundary 
proceed in opposite directions on the X plane. The behavior near the real 
root boundary involves four roots, two of which meet at 1/rev and two at 
-1/rev; hence, the roots occur as complex conjugate pairs always. 

In forward flight (y > 0) the behavior of the roots for individual blades 
in the rotating frame is the same as that of the hover roots (i.e., as shown 
in fig. 9(a)), just with an order y 2 shift. Since the blades are independ- 
ent, the transformation to the nonrotating frame can only shift the locus by 
ii in forward flight as it did in hover; so figure 9(b) presents the behavior 
of the roots of individual blades in the nonrotating frame for forward flight 
as well as for hover. The three-bladed gimballed rotor, however, exhibits 
different behavior near the real root boundary in forward flight. Instead of 
the two roots meeting and then proceeding in opposite directions at 1/rev fre- 
quency, the two branches of the loci only pass close to each other as illus- 
trated in figure 9(c). That is, as y increases, the locus intersection 

pulls apart. In terms of the y locus at a given y > 0, as the branches 

approach 1/rev there is a ± order y 2 change in the damping (with a corre- 
sponding change in the conjugate roots, so this is not critical region 

behavior); eventually the roots transition to more like the two real root 
behavior (at 1/rev frequency nonrotating), but with a ± order y 2 differ- 
ence in the frequency. Such behavior is, in fact, typical of the dynamics of 
coupled degrees of freedom; the root loci do not cross, but rather only pass 
close to one another, pulling farther apart as the coupling increases (in this 
case, as y increases). The source of this behavior is the last term in the 
radicand of equatiofi (104) for X. There is, in addition, the order y 2 
change in the ±i term of equation (104), so all this occurs really at not 
quite ±l/rev. 

It follows then, that when y > 0 the three-bladed gimballed rotor does 
not have a definite real root boundary. Rather there is a gradual transition 
from primarily complex root behavior to primarily real root behavior. Examine 
now the locus in the neighborhood of the r eal root boundary. Write 
y = Y 0 + y 2 y where y Q = Kp + /v 2 + Kp 2 , so y is order y 2 from the hover 
real root boundary. Look at the root near the 1/rev frequency, and near the 
-y/16 real part; that is, let AX be defined by 

X - ijl ♦ u 2 £ [j- Qj) ♦ v“JJ - ^ ♦ AX (106) 

Then, from equation (104) , 
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This may be separated into real and imaginary parts: 


(ReAA) 2 - (ImAA) 2 


y2 [x(lS " k p) ‘ T Kp ‘ T v<t ] C108) 


(ReAA) (ImAA) = - |y 2 ^v 2 


(109) 


The last equation shows that the locus on the AA plane (i.e., on the A 
plane, but with a shifted origin) is a hyperbola. The point of closest 
approach of the two branches occurs at 


ReAA * -ImAA = ±]x 



( 110 ) 


and the minimum separation of the branches is there 


2 1 AA | 


mxn 


8 Ao 

M l/l6 v 


(HD 


The separation of the branches then is order y. This point of closest 
approach may be taken as the definition of the real root boundary for this 
case. With ReAA = -ImAA, equation (108) then gives the boundary as 


or 


Y* 


IS 



Y 0 „ 
lS K P 



2 f v2 (~ K P * * K P 2 ) * K p 

/v 2 ♦ Kp 2 


( 112 ) 


For Kp * 0 this reduces to y/16 * v[l + y 2 (8/9)u 2 J. Notice that this is the 
same as would be obtained if the last term j • Uie radicand of A (eq. (104)) - 
which is causing all this behavior - were s^:> < ’ ignored. Setting the 
radicand to zero without this term gives exa ’ the above boundary. The 
boundary obtained is similar in behavior to that of the individual blades; but 
again the second order u 2 in the radicand of equation (104) has only half 
the magnitude. Hence, the Ay of the boundary for Kp * 0, for example, is 
only half as large as for the independent blade. 
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Near Vrev Frequency 


Consider now the case ImXg = 1/2, when the nonrotating roots are at 
3/2/rev and -Vrev frequency. The order y secular equation for Bqi- is 
unchanged for this case, hence, the low frequency mode (Xq " 1) does not have 
a critical region when the rotating frequency is near H/rev. The order y 
root is the same as equation (98) then, that is, there is no forward flight 
influence on the lower frequency root to order y when ImXo * 1/2. The 
secular equation for Boi+ i s » however, changed for this case, and so the 
high frequency mode does encounter the Vrev critical region. The complete 
rotor behavior is given, of course, by both the low and high frequency modes 
taken together; when the high frequency roots encounter the Vrev region, the 
entire rotor does, even though only two of the four roots participate in the 
critical region behavior. 

When ImXo * 1/2, the order y secular term for $oi+ is 



(113) 


(114) 


This is the same secular equation and D 2 as for the Vrev critical region of 
the independent blade. It follows then that the root behavior is the same as 
that solution, specifically the expressions for the critical region boundary, 
and for X near the region (eq. (32) and (34)) are applicable here as well - 
adding i to X to transform it to the high frequency mode in the nonrotating 
frame. 


Near 1/rev Frequency 

Consider the case ImXo ■ 1, so the nonrotating roots are a 0 and 2/ rev. 
The order y analysis is applicable, since ImXo t 1/2, and the order y 2 
analysis is correct for the high frequency mode even when ImXo * 1. Hence, 
the high frequency roots (near 2/rev frequency in the nonrotating frame) are 
given by the same expression as for the roots away from the critical regions 
(eq. (104)). The low frequency mode, however, (near 0/rev nonrotating) has an 
additional term in the order y 2 secular term when ImXo * 1* instead of 
equation (101), Bn- is now given by 
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)-^r + t (x ° + Kp)6u -] 


s |( 2X ° + t) 
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3<l»2 




Y 

(X 0 + ^p) + 


♦ $■ (» 0 ♦ 2Kp)(-i)A.]6 02 -)e X ‘*‘ 

♦ (1 ♦ Kp)9o 2 -]e Tl * 1 < U5 > 


Unless Ai is real, so Ai * A lf the secular equation for 6 0 2-O2) is the 
same as before, however. Then, equation (104) holds for the lower frequency 
mode also, except when the frequency is order y 2 from 1/rev (rotating). 
Requiring Ai be real means y\ = 0; the secular term of equation (115) is 

then 



which has similarities to the results for both the individual blade and the 
teetering rotor. The boundary of the critical region is defined by D * 0, 
letting y * yq + Ay with Ay order y 2 small, then the boundary is 


Kp) ■ k 2 (C, ± C 2 ) (118) 


or in terms of the comer u 


48 


oy no „ \ 
2 ..-2 - ig ^ 


Wl ,W2 


Cj iCj 


( 119 ) 


This is the usual form, with now 


r I /lif /Ii\ 

Ll = " 2 \12/ ~ \ 6 } 

C 2 = A * Kp 2 


2v2 'f "p 


♦ 4K p 2 
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if Y ° 

K P T5 


( 120 ) 

( 121 ) 


C 2 , which gives the width of the critical region on the y ~ v* plane, is the 
same as for the teetering rotor; hence, the width is, in general, smaller than 
for the individual blade, except for v « 1, in which case the widths are 
equal. Cj , which gives the offset of the region, is similar to the individual 
blade result. The first term is new, however, corresponding to the new order 
p 2 shift in the ±i frequency of the roots away from the critical region 
(eq. (104) j; and the second term has half the magnitude, corresponding to the 
second order u 2 tern in the radlcand of equation (104) . The eigenvalue of 
the low frequency mode is now 


• - & t • " 2 1 ^ • l) ] 1 1 U (t? • 2lt p) /( l ' T^X 1 ' ~§) 

( 122 ) 

So the eigenvalue also has a new term, an order u 2 damping change, 
corresponding to the third order y 2 term in the radicand of equatio n (104) . 
For Kp - 0, the criterion for ImXo = 1 reduces to y 0 /16 = Ak - 1, 
equations (120) and (121) become 

Cl * ‘ 7 (if) ' (t) h “ ‘ 77 (v2 • 1)(v2 * 3) (123} 

Y j 

C 2 ■ - J /j 2 - 1 (124) 

so the boundary reduces to 
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( 125 ) 
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and the eigenvalues to 

x e - T5 t 1 ' y2 I (v2 ‘ x) ] 



(126) 


Compared to the individual blade result, the boundaries are identical if 
v * 1, but if v > 1 this boundary is offset more, with a smaller width. Com- 
pared to the teetering rotor, the width is same but the teetering rotor is 
not offset even when v > 1. 


Summary 

For the hover limit, y * 0, the four eigenvalues of the three-bladed 
gimballed rotor are 

X = ' TS + i /> 2 + } K P - (t?) (127) 

and their conjugates. For the hover frequency away from %/ rev or 1/rev 
(rotating), the roots to order y are the same as the hover roots, and to 
order y 2 , they are given by equation (104) above. The influence of forward 
flight consists of: an order y 2 frequency change, similar to the individual 

blade but with reduced magnitude; an order y 2 change in the ±i frequency 
shift due to the transformation to the nonrotating frame; and a ± order y 2 
change in the damping of complex roots of the frequency for real roots. The 
latter change has the effect of pulling apart the real root intersection of 
the locus (at ±l/rev frequency nonrotating) , so that the root loci only pass 
close when y > 0 rather than actually intersecting (the separation is order 
y) . The criterion for divergence stability is, in this case 



The tyrev critical region is identical to that of the individual blades. 
Only the high frequency mode, at 3/2/rev nonrotating frequency, participates 
in this region, however. The low frequency mode, at ij/rev, has the same 
behavior as it does away from the critical region - that is, no forward flight 
influence at all to order y. 


50 



If Ay * y - Y 0 is order u 2 sm all, where y p is the value required 

for ImXo 3 1, that is, y 0 /16 = Kp + /> 2 + Kp 2 - 1; then the low frequency 
mode, at 0/rev nonrotating frequency, encounters the 1/rev critical region. 
The high frequency mode does not participate in the critical region behavior; 
its roots are given by the same expression as away from the critical region, 
that is, equation (104). The low frequency roo«.j> are given by 


= • lk [ l - p2 1 (v2 ■ 1} ] 1 1 l£ G? - 2 k p)/( 1 ■ ■ fi) 


where the comer u 2 "ire 
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with Ci and C 2 given by equations (120) and (121) above. The behavior of 
the loci in the critical region is similar to what has been seen before with 
the individual blades and the teetering rotor. There are now, however, order 
u 2 changes in the frequency and damping (in Ci, and explicitly in ReX 
above) corresponding to the changes seen away from the critical region, as 
well as some changes in the magnitude of the forward flight influence. 

The y - w plane boundary for the two real root region (il/rev 
intersection nonrotating), writing y * Y 0 ♦ Ay, I s 


Ay 

ir 
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2 1*1 

16 


v 2 (-Kp ♦ /v 2 + Kp 2 ) ♦ Kp 


(131) 


where y 0 /16 » K p ♦ /v 2 ♦ Kp 2 . This is only a soft boundary, however, the 
point of closest approach of the branches when u > 0. The '-s/rev critical 
region boundaries (3/2/rev nonrotating) are 
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Ay . IT 




Kp ♦ 4Kp 4 


Y 0 

IS 


(132) 


- K, 


where y 0 /16 » Kp ♦ /v 2 ♦ Kp 2 - (1/4); and the 1/rev critical region 
boundaries (0/rev nonrotating) are 




( 133 ) 
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where y 0 /16 = Kp ♦ A> 2 + Kp 2 - 1. For Kp » 0 these boundaries reduce to 


- u 2 - v 3 
16 y 9 


with rr ■ y 
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respectively. 


GIMBALLED ROTOR, FOUR BLADES 


The equation of motion for a four-bladed gimballed rotor are 



The only forward flight effects remaining are order y 2 , the order y moment 
cancelling internally at the hub as they did for the teetering rotor. The 
only periodic coefficients remaining are order y 2 , hence, only a 1/rev 
critical region is expected. 

The hover result i:» identical to that of the N ■ 3 case; that result, 
in fact., holds for all N >. 3. Only order y 2 terms appear in the equations 
of motion, so an expansion in y 2 is used: 
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_L „ _J_ + y 2 _L- + . . 

34) 34»0 ^2 


g = 6 0 + U Z &2 + 


v * y 0 + v y 2 + 


• • • 


The order 1 solution is the hover limit again, so is the same as the 
result. The solution for 80 is 


I. - Re[ 6 „ 2 .C* 2 ,...)Qe (X »* i) *» * 6o 2 .t* 2> . • •)(-;>—] 


(138) 


Order u 2 Results 
The order y 2 equation is 
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Assuming that ImXg j 1, the secular term is 

30O2± 


d<l>2 


- X 2602 ± * 0 


(139) 
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where 

y 2 Yn 

X (*0 ♦ Kp) + -f Kp 

* 2 " 2iImX 0 

The solution is 6o2± * &04±(^4»* • «)e* 2< ^ 2 . Hence, the eigenvalues to order 
y 2 are 


X * Xq + U 2 X2 ± i 

m tj - 15+ i / 2 + ^ + y2 > 8 K P ‘ ^ 141) 

and the conjugates. This is the same result as for the teetering rotor roots 
away from the 1/rev critical region. It follows then that the real root 
boundary (at ±l/rev nonrotating) and the divergence stability criterion are 
the same as for the teetering rotor. 


Near 1/rev Frequency 

If the hover roots are order y 2 from 1/rev frequency (rotating) - that 
is, ImXo * 1 - then, the periodic coefficients contribute to the order y 2 
secular terms: 
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36 n?+ . Yn - 

~9^ ^2802+ ' 1 32 ^ + K P^02± = 0 

Comparing this with the teetering rotor 1/rev secular equation, it follows 
that the four-bladed gimballed rotor has identical expressions for the roots 
and boundaries of the 1/rev critical region. 


Summary 

The behavior of the four-bladed gimballed rotor is the same as that of 
the teetering rotor flap stability - except, of course, that the gimballed 
rotor has four roots, in the nonrotating frame, hence, shifted by ±i from the 
teetering rotor roots (rotating). As for the teetering rotor, it also follows 
here that the constant coefficient approximation to the equations of motion, 
even in the rotating frame so the only forward flight influence retained is 
the factor of (1 + y 2 ) in the mean of Mg, gives exactly the correct 
eigenvalues except near the 1/rev critical region. 

Away from the critical region, the eigenvalues of the four-bladed 
gimballed rotor are 


A - ±i - + i y4 2 + (1 + y 2 ) | K p - (143) 

and their conjugates. There is no Vrev critical region. The 1/rev critical 
region is the same as for the teetering rotor; both the high frequency and low 
frequency modes (at 2/rev and 0/rev nonrotating) participate in the critical 
region behavior. The y - y plane boundaries are the same as given for the 
teetering rotor. 


GIMBALLED ROTOR, FIVE OR MORE BLADES 


For a gimballed rotor with five or more blades, the equation of motion is 





( 144 ) 
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So to order y 2 there are no periodic coefficients when N 2 . 5. All the 
periodic flap moments cancel internally at the rotor hub, and do not contrib- 
ute to the net pitch and roll moments on the rotor disk. The result is a con- 
stant coefficient differential equation for the gimballed rotor with five or 
more blades. The eigenvalues of this rotor are given by the roots of the 
characteristic equation of 


2X - * (l - $) 


2A - I 


X 2 + I 
A 8 



or 

[a 2 + A + v 2 - 1 + Kp ^ (1 + y 2 )J + (2A + ■g) ~ C 1 + K P 2 ) * 0 

(145) 

As a constant coefficient equation, the eigenvalues are simply the roots of a 
polynomial, although in this case it is a fourth-order polynomial, which may, 
in general, only be solved numerically. Since the roots for the limit y ■ 0 
are known, however, - that is, the hover roots - a perturbation technique may 
be used to obtain explicitly expressions for A including the influence of 
forward flight. 


Only order y 2 and y 4 terms appear in the characteristic equation, so an 
expansion in y 2 is used. Write A * Aq ± i + y 2 A 2 + . . . and the conju- 
gates, where Aq is the rotating hover root as usual: 

x o a + i / v2 + T k p ‘ (l&) (146) 

The Lock number is also expanded as a series in y: y a Yq + U 2 Y 2 * • •• The 

order 1 characteristic equation then has just the solutions A = Aq ± i, 

Xq ± i. That is, the order 1 term in the expansion of A above is indeed the 
correct hover limit. 


Order y 2 Results 

The order y 2 terms of the characteristic equation are 

F( x 0 - i) 2 + ~j^ ( x 0 * i) + v 2 - 1 + Kp "^J^CAo ± i)A 2 A 2 

♦ ^ (A 0 ±i) + Kp ^■♦Kp ^]+|(A 0 ±i) *^](2X 2 +x)" 0 


56 



or 


(ImX 0 


± l)[(2i 


ilraXg) X 2 + Kp 


1L 

8 




(147) 


Assuming that ImXg / 1 when the lower frequency root Xg - i is being 
considered, so ImXg -1/0, then the solution is 


X 2 = 


(x 0 + Kp) + y -f Kp 
2ImXg 


(148) 


It follows then that the eigenvalues to order p 2 are 


X 


* ±i - 




(149) 


and the conjugates. The real root boundary and the divergence criterion are 
the same as given for the teetering rotor. 

Notice that this is the same result as would have been obtained if a 
constant coefficient approximation were made in the rotating frame before 
finding the net pitch and roll moments. Actually, the expansion of X as a 
series in p 2 is not really necessary, for to order p 2 the last term in the 
characteristic equation (the u 4 term in eq. (145)) drops and X may be 
found directly, exactly as for the hover case. However, such a procedure does 
not show that the solution is not valid for the low frequency mode near 1/rev 
(rotating frequency) . 


Near 1/rev Frequency 

This is where the critical region is encountered when the equations have 
periodic coefficients, that is, when ImXo = 1. The high frequency mode is 
at 2/rev nonrotating frequency, and the roots must remain complex conjugates 
since this is a constant coefficient equation. Hence, the 1/rev critical 
region behavior cannot be encountered for this mode. Indeed, the order p 2 
result obtained above holds for the Xo + i root, even when ImXg = 1 (the 
upper sign is used in eq. (147), so the first factor has the value 2). The 
low frequency mode, however, is at 0/rev nonrotating frequency, so these roots 
are able to meet and proceed in opposite directions along the real axis. This 
is the 1/rev critical region behavior which is allowed for this constant 
coefficient equation because the transformation to the nonrotating frame puts 
these roots near the real axis instead of 1/rev. 

Consider the low frequency mode when ImX 0 = 1; that is, 

X a -(Yo/ 16) + p 2 X 2 + . . .. The order p 2 characteristic equation is 
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identially zero then (eq. (147), using the lower sign so the first factor has 
the value 0 when ImX 0 = 1), and it is necessary to go to order u 4 to find 
* 2 . The order u 4 terms in the characteristic equation are for this case 
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Setting D 2 = 0 gives X 2 = so X *» -y/16 to order u 2 . Hence, D 2 * 0 

is the boundary for two real roots in the nonrotating frame, or the 1/rev 
boundary in the rotating frame. In general, the root is 


X =* 



■-**»*» 


(151) 


The expressions for X and D are identical to the corresponding ones for the 
teetering rotor, or for the four-bladed gimballed rotor, near the 1/rev criti- 
cal region. Thus, the gimballed rotor with five or more blades encounters a 
1/rev critical region, although only the lower frequency mode (near the real 
axis nonrotating) participates in the critical region behavior. TV ' boundary 
of the region, and the solution for X, are the same as were given or the 
teetering rotor (-i to get to the low frequency nonrotating mode) . 


Summary 

The gimballed rotor with five or more blades is described by constant 
coefficient differential equations. The behavior of the roots is, however, 
nearly identical to that of a teetering rotor or a four-bladed gimballed rotor; 
the latter is perhaps a better comparison since it is also a case with non- 
rotating degrees of freedom. Away from 1/rev rotating frequency, the roots 
are given by 

X = ±i - fa ♦ i/v 2 ♦ (1 ” u 2 ) $ K P - (jfcf (152) 
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and the conjugates. The real root boundary (±l/rev nonrotating) and the 
divergence criterion are the same as for the teetering rotor. 

A critical region is encountered by this rotor, when the rotating 
frequency is near 1/rev. The expressions for the boundary and the roots near 
the critical region are the same as for the teetering rotor. Only the lower 
frequency mode participates in the 1/rev critical region, however, for it is 
the fact that this mode is near the real axis (nonrotating) that allows the 
critical region behavior to occur with a system described by constant coeffi- 
cient differential equations. The high frequency modes are near 2/rev for 
this case, and so are still given by equation (152), that is, have the same 
behavior as away from the critical region. 

The equations of motion for the gimballed rotor with N >. 5 are identical 
to the constant coefficient approximation of the equations for the three or 
four-bladed gimballed rotors (i.e., eq. (87) or (137) with the periodic terms 
dropped) . Hence, the solution for N >. 5 may be considered as a constant 
coefficient approximation to the dynamics of the N = 3 or 4 cases. 

As a constant coefficient approximation to the N = 4 case, this present 
solution gives exactly the correct roots with the exception of the high fre- 
quency mode near the 1/rev critical region. These roots for the N * 4 case 

encounter the critical region along with the lower frequency roots; but for 

the constant coefficient approximation s-. -',h behavior is not allowed for the 
high frequency roots if they are to remain complex conjugates. The behavior 

of the low frequency roots is given correctly everywhere, including the 1/rev 

critical region, because in the nonrotating frame that behavior occurs on the 
real axis. Away from the 1/rev critical region, it is, in fact, possible to 
make the constant coefficient approximation in the rotating frame, before 
finding the net pitch and roll moments, and still obtain the correct 
expressions for the roots of the four-bladed gimballed rotor. 

As a constant coefficient approximation to the N = 3 case, this present 
solution is not really good anywhere. The errors are not even order y 2 only, 
because the three-bladed gimballed rotor encounters the Vrev critical region 
where there are order y effects of forward flight. The order y 2 effects 
on the roots, both away from and near the 1/rev critical region, are, in addi- 
tion, quite different for the N = 3 case. 


EQUATIONS OF MOTION IN THE NONROTATING FRAME 


Consider a rotor with N independent blades, each with rotating natural 
frequency v. The flap motion of the entire rotor is described by the N 
degrees of freedom 8( m ), and the N rotating equations of motion given above 
(eq. (57)), in contrast to the gimballed rotor, which is described by onl / two 
degrees of freedom and equations for all N _> 3. As for the gimballed rotor, 
however, the motion of the rotor with independently flapping blades may also 
be described in rhe nonrotating frame. The following new degrees of freedom 
are introduced: 
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^ ' 


where 4^ = + mAtp is the azimuth location of the mth blade and 

A 1 J 1 = 2n/N is the interval between the blades. The flap motion of the mth 

blade is then given by 


gW = Bo + £(&nc cos "'I’m + &ns sin "I’m) + e N/2 t -1 ^ m ( 154 ) 

where the sum over n goes from 1 to (N - l)/2 for N odd, and from 1 to 
(N - 2)/ 2 for N even. The & N / 2 degree of freedom only appears if N is 
even . 

This is a Fourier coordinate transformation from the N degrees of 
freedom describing the rotor motion in the rotating frame to the N 
degrees of freedom Bo» 3 nc , 3ns > &N/2 describing the rotor motion in the 
nonrotating frame (for further discussion, see, i.e., ref. 13). The Bo 
variable is the rotor coning motion; 3ic and $is are the tip path plane longi- 
tudinal and lateral tilt degrees of freedom as for the gimballed rotor. The 
Bn /2 motion is similar to the coning except that the blades alternate in up 
and down motion. This coordinate transformation must be accompanied by a con- 
version of the equations of motion from the rotating frame (eq. (57)) to the 
nonrotating frame. This is accomplished by operating on the equations of 
motion with the following summation operators: 

i £(•••). §£(...) cos n* m , £ E(...)sin wlfc, £ 
mm m m 


(155) 


The summations combine the equations for flap moment equilibrium into the 
rotor moment equilibrium appropriate to the nonrotating degrees of freedom. 

For example, for the coning motion, the first onerator finds the coning moment; 
and for the Bic mid Bis motion, the operators fine che net pitch and roll 
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moments on the disk as was done for the gimballed rotor. The resulting N 
equations in the nonrotating frame are coupled, even for the case of a shaft 
fixed rotor as here in contrast to the rotating equations, which are not 
coupled at all (eq. (57)). 


The result of the summation operators depends on the number of blades (N) 
where, as here, the equations have periodic coefficients. Hence, the set of 
N nonrotating equations depend on N. The complete sets of differential 
equations describing the rotor flap motion of N independent blades excited 
by blade pitch inputs only are given below for the cases N = 3 and N = 4. 
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The derivation and discussions of these equations may be found in the 
literature, for example in references 13 and 18. 

The transformation to nonrotating degrees of freedom and equations of 
motion has the effect of sweeping the periodic coefficients from the lower 
degrees of freedom, especially as N increases. There are, however, always 
periodic coefficients in the nonrotating equations where y > 0, just as there 
are in the rotating equations, for the same system is being described and a 
constant coefficient set of equations could never give the behavior that has 
been found with the periodic coefficients. 

The flap dynamics described by these sets of differential equations is 
nothing new. This is still a set of N independent blades and simply observ- 
ing the motion in the nonrotating frame does not change the nature cf the sys- 
tem. Hence, the eigenvalues must still be as calculated above for the case of 
a single independent blade; the only change is a ±ni for the 8 n _/6 ns modrs, . 
that is, a ±n/rev shift in the frequency to account for the transformation 
from rotating to nonrotating frames. The perturbation solution for these 
equations, including the influence of the periodic coefficients in forward 
flight, has already been obtained then. The value of the description of the 
rotor by these nonrota';ing degrees of freedom and equations lies in its use for 
dynamics involving the helicopter body or shaft motion, aerodynamic gust, or 
any other excitation from the nonrotating frame. The rotor responds to such 
excitation as a whole, in nonrotating modes of motion, so this description of 
the rotor motion and moment equilibrium is appropriate for studying such prob- 
lems. For the current investigation, the use of these nonrotating equations 
is that they provide a basis for a constant coefficient approximation for the 
rotor flap dynamics. 


CONSTANT COEFFICIENT APPROXIMATION 


As should be quite apparent by now, differential equations with periodic 
coefficients require considerably more analysis than constant coefficient 
equations, even to simply find the eigenvalues. Constant coefficient differ- 
ential equations are much preferable for the studies of dynamic systems. Can 
the fli dynamics of a helicopter rotor in forward flight be adequately 
describeu by some constant coefficient approximation to the equations of 
motion? The constant coefficient approximations for the teetering and gim- 
balled rotors were discussed along with their perturbation solutions including 
the periodic coefficients; this section will be concerned with the case of 
individually flapping blades. 

Consider the rotating equation of motion for 0( m ) (eq. (57)). The 
constant coefficient approximation to that equation, that is, using only the 
mean values of the coefficients, is 

0 ♦ J 0 ♦ [v 1 ♦ (1 ♦ u 2 ) J Kp]e = 0 (158) 
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The eigenvalues of which are 


1 ■ - TS ! i/* 2 * » * ^ l IV - (n) 2 (159) 

The only influence of forward flight is the order y 2 increase in the mean of 
Mq, so in the effectiveness of Kp. Consider this equation and its roots as a 
representation of the flap dynamics of an independent blade in forward flight. 

A correct estimate of the eigenvalues to order y is obtained with such a 
representation, except near the h/rev critical region, simply because neither 
this constant coefficient approximation nor the correct solution away from the 
Vrev region show any order y effects. To order y z , this approximation is 
not correct even away from the critical regions. 

With the constant coefficient approximation made in the rotating equation, 
the result is virtually just the hover root; so the question of the applica- 
bility of that approximation is really whether the hover roots may be used as 
a measure of the flap dynamics in forward flight as well as in hover. This 
approximation misses the Vrev and 1/rev critical regions entirely, and so 
provides no information at all about those effects. Away from the critical 
regions, however, the correct roots show only quite small influence of forward 
flight, as for example case (c) in figure 1. This constant coefficient 
approximation - that is, the hover root - may be considered a reasonable 
representation for such cases. It should be noticed, however, that the 
influence of the critical regions becomes greater as y increases, so this 
approximation ir. st eventually break down for every case. 

The nonrotating equations of motion provide another source of a constant 
coefficient approximation to the flap dynamics of a rotor with N independent 
blades. The transformation to the nonrotating frame is accomplished first, 
and then the averaged coefficients are found. The harmonics of the coeffici- 
ents of the rotating equation contribute to the constant coefficients in the 
nonrotating frame, hence, this procedure retains more of the influence of the 
periodic coefficients. For N * 3, the resulting equations of motion are 


64 



Bic I + 0 


0 Y 2C 

-yM. 4 m. 

B 2 B 


- 4 m 


-yM. rr M. 


1m! c 4h! ! 


+ -yM c V* - 1 


4m? s 


V 2C 0 V 2( 

- 2 M e s 


yM? - 1 M? C 
8 1 B 


0 Y 2C 

yM 6 ♦ i m 9 


®0 " Kp®o 
01C - KpBic 
pis - K P 6 1Si 


( 160 ) 


where use has been made of the following Fourier series representations of the 
flap moments in forward flight: 

0 is 2c 

M 0 = M 0 + Mg sin ^ ♦ Mg cos 2i|> ♦ . . . 

0 is 2C 

M. » M. ♦ M. sin i|> + M. cos 2i|> + . . , 

B a 6 6 

1C 2S 

Mg « Mg cos ij» ♦ Mg sin 2 i|i ♦ . . . 


The missing harmonics is a general result valid for all y, even including the 
reverse flow region effects. Substituting for che harmonics of the flap 
moments from equation (2), that is, neglecting reveits flow as usual for the 
present order y 2 analysis, the equations are 


•* .4- 




(161) 


This is the constant coefficient approximation, in the nonrotating frame, for 
the flap dynamics of a rotor with three independent blades. It may also be 
obtained from equation (156), by dropping the 3/rev terms in the coefficients. 

The approximate eigenvalues are then solutions of th • characteristic 
equation 

*1 2 

^> 2 + j X ♦ v 2 + Kp (1 ♦ w 2 j]||x 2 ♦jA + v 2 -l+Kp|-(l + u 2 )^ 

* ( 21 *i) 2 - ""(li s) 2(I * i 'p 2) | • 1,2 (tt)' 2 '* * 2K p>[ 2x *i 0 -t)] 

- U 2 (^) 2 2(».2Kp)2[x2.J » - 1 J (l .0 (1*2) 

Only u 2 and u 1 * terms appear, so consider an expansion in u 2 : 

X * \ 0 ♦ , 2 x 2 ♦ . . . , or X » X 0 ± i ♦ U 2 X 2 ♦ . . ..and the conjugates, for 
the coning mode and Bic/Bis modes, respectively. The Lock number is also 
expanded as a series: y * yo * U 2 Y 2 4 • • •» and Xq is the rotating hover 

root as usual: 

x o * • if ♦ i/v 2 4 x k p ' 

To order 1, the solution for the roots is just X » Xq, X 0 ± j and their 
conjugates, verifying that the above expansion is correct for the hover limit. 


Coning Modes 

Consider the coning mode roots in forward flight, that is, the expansion 
X » Xq + y 2 X 2 ♦ . . . • The order u 2 terms in the characteristic equation 
give then 
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(>x 0 x 2 * ^ X 2 * T ( 2X ° * $ 1 

.(!0) ? 2(X 0 *2Kp)(2X 0 *^)*(f“) 2 2(X ll *2 l Cp) 


2 _ 


= 0 


[1 - (2imX 0 ) 2 ][2iI°H(‘2*lf) * T (- TS* k p) + T Kp ] 

* (S)H 2 - T K P * n * 2 ) - 0 

Assuming that IuX„ t 1/2, the root to order u 2 is then 


(163) 


X = 


- ± i / 2 + (1 + V* 2 ) g K P ‘ (jk) 


t 2 H 

1 + P 9 


v 2 -lK P + 4K p 2 


v2 + 8 K P 




(164) 

™ - «,* h LSSS t ^r t ^ r^t^rTaL^e^r^l 

foi this root are also correct them This “/ifmissefthe 

the coning mode approximate root when ImX 0 - 1/rev, however, 

1/rev critical region. 

Now consider the 

^rcha^i^'eUion aS 

the order u 2 term is 


(2X0 *,t ^ X, X„tKp ^[-Z^XoXi M *T Xo^P 

* 2 ( 2 X 0 * x)( 2X ‘ *t)] " (it) 2(^0 * 2Kp)(x 0 * T 


- 2Kp 


■)■ 


or 


(e. * fe)’ ■ -[^(ft - >r)]’ • (S/(- - ? 't • 


= -D 2 


(165) 
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The root to order y is then 


X = Xq yXj = - 




( 166 ) 


These are exactly the same expressions for D and X as were found for the 
%/rev critical region by the perturbation analysis in the rotating frame, 
including the periodic coefficients. So the coning mode of this constant 
coefficient approximation encounters the %/rev critical region, with the 
correct boundary and roots. Now the roots in a critical region have a plus 
and minus increment in the damping due to the periodic coefficients, that is, 
the roots are not complex conjugates. What is happening to allow the roots of 
a constant coefficient approximation to exhibit this critical region behavior 
is the following. The low frequency mode roots Xq - i and the conjugate are 
also at ±%/rev when IjoXq = 1/2. Hence, with the coning and low frequency 
mode roots there are a total of four roots, two at %/rev and two at -%/rev, 
which must participate in this %/rev critical region behavior. Hence, the 
critical region behavior can occur while the roots remain as complex conjugate 
pairs as is required for constant coefficient equations. 

High and Low Frequency Modes 

Consider the low frequency and high frequency rotor modes, that is, the 
expansion X = X 0 ± i + y 2 X 2 * • • •• Then the order y 2 equation is 

C ImX o ± 1) j(+2ImX 0 - 1) [c2iImXo)X 2 + (> 0 + Kp) + Kpl 


■(»)' 


(X 0 t i + 2Kp)(X 0 + 2Kp)| = 0 (167) 


Assume for now that ImXo t 1/2 or 1 for the low frequency rotor mode; then 
The roots to order y 2 are then 




(169) 


and the conjugates. This result is similar to what was found for the three- 
bladed gimballed rotor away from critical regions, so the behavior is familiar 
from that discussion. The only change, in fact, is a reversal of the signs of 
the ± order y 2 terms, hence, nil the analysis of that case may be used 
here. 


The main feature of the behavior is that when y > 0 the real root 
intersection of the y loci - at ±l/rev nonrotating - pulls apart. The order 
y 2 changes in the frequency and damping are in the opposite direction from 
what is shown in figure 9(c). The closest approach of the two branches has an 
order y width: 

2 Mmin - ^ f /}% v (170) 


where yq/ 16 = Kp + /v 2 + Kp 2 is the y for the boundary when y * 0. 
Taking the point of closest approach as the real root boundary, it is then 



^ v 2 (-Kp + /v 2 + Kp 2 ) + Kp 
/v 2 + Kp 2 


(171) 


The criterion for divergence stability for these four roots is 


(1 + u 2 ) | Kp > -v 2 


1 + 



(172) 


There is this important difference from the case of the three-bladed 
gimballed rotor: this solution does not represent the behavior of any real 

system, but rather is being considered as an approximation to the rotor flap 
dynamics in forward flight. As such, it has then order y 2 errors in both 
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the frequency and damping compared to the correct solution (eq. (24)). Away 
from critical regions, the correct solution shows only a small influence of 
forward flight anyway, so these errors are probably not large either. 

This behavior of the roots n.ar the real root boundary has been seen 
twice now: with the three-bladed gimballed rotor, and here for the constant 

coefficient approximation to the independent blade dynamics. The origin of 
this behavior is the order y, 1/rev term in the flap damping M. . That term 

gives the 3/rev periodic coefficients for the gimballed rotor which produce 
the term in the radicand of X that gives this behavior. For the constant 

is 

coefficient approximation, the flap damping harmonic M. gives terms cou- 

6 

pling the 6i c and gjs motion, producing the observed behavior. With inde- 
pendent blades observed in the rotating frame, the loci intersect at the real 
axis as is characteristic of root loci; observed in the nonrotating frame, the 
behavior is shifted by ±i so it occurs at ±l/rev now, but the blades are 
still independent so the root loci still intersect. The root loci of coupled 
degrees of freedom do not show that behavior, however; rather they simply pass 
close to one another, the minimum separation being a measure of the magnitude 
of the coupling. The gimballed rotor couples the blades by requiring that the 
rotor move as a whole, in 6i c and 3 ls motion only. The constant coefficient 
approximation couples the blades by dropping terms from the equations of 
motion describing the independent blades, creating a description of some new 
system with coupled degrees of freedom. (How closely the dynamics of this new 
system might represent the rotor with independent blades is what is being 
examined here.) 

When ImXg = 1/2, that is, when the hover rotating frequency is near 
Vrev, the low frequency rotor mode (at ~Vrev non. '"’•ting frequency) 
encounters the Vrev critical region. It is not nt sary to go looking for 

this solution, because the only way the coning mode uld exhibit the critical 

region behavior with a constant coefficient equation is if the low frequency 
mode joins it, so the roots may remain complex conjugates even inside the 
critical region. The critical region behavior found is the same as that of 
the perturbation solution including the periodic coefficients. So the con- 
stant coefficient approximation gives the correct behavior of these four roots 
near the Vrev critical region. The high frequency rotor mode (at 3/2/rev 
nonrotating frequency) does not, however, participate in this behavior. 

Indeed, the above result (eq. (169)) is still valid for the Xq + i root even 

when ImXo = 1/2; the error in that root is then order y. 

Now consider the low frequency mode (Xq - i) when ImXo = 1; these roots 

are at 0/rev nonrotating frequency, that is, to order 1 X * -y 0 /16. The high 

frequency mode gives no problems in the order y 2 characteristic equation 
even when ImX 0 = 1, so equation (169) gi/es the roots there. For the low 
frequency mode, however, the order y 2 characteristic equation (eq. (167)) is 
identically zero when ImX 0 = 1, so it is necessary to go to order y 1 * to 
find X 2 . The order y 1 * characteristic equation is, when ImXo = 1: 
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or 

»« * IS * - 7 fef( jt ' 2k p) 1 10 

where 


D 2 = 

/Y„ \ y 9 /Yo\ 2 v2 '^ Kp + 4Kp2 1 / Y o\ 2 v. Y ° 

p- K ») w* (t) 8 -2 (n) - Kp TsJ 

(Sfjt *T (S- 2Kp )] 2 *[ Kp - IT (S- 2Kp ) 2 ] | 


and the root to order y 2 is 


- * [' * » 2 if (il 2ltp )] 4 11,20 


This is the critical region behavior again, possible with the constant 
cient equations because for this mode the nonrotating roots are at the 
The critical region boundary is given by D =0, or the corner 


axis. 




fro „ \ Ay 
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where here 


Ci = - 



v 2 - x k p + 4K P 2 
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c * ’ ft j[ *t(t6- 2k p)] 2 4 [ K P - ft (if - 2l< p) 2 ] 2 


'ft /[' 4 T ^ ' ‘I]' 4 [ K P * T t'' 2 - »(ft - 2K pjf 


(178) 

and yq/ 16 = Kp + /v 2 + K p 2 - 1 for the 1/rev region as usual. The root may 
be written then 


> - £ [l * u 2 f (V 2 - 1)] ± i fHfe - Kp)/|l 

■ £)( ' i?) 

(179) 

and the boundary on the y - V plane is 



fiy /l", 

16 W Y ° K 

16 " K P 


(180) 

For Kp - 0, these reduce to Y 0 /16 = /v 2 - 1, 



c, - - | (v* - 1)2 


(181) 

n ^ 0/1 32 r 0 n , 256 2 r 2 

C2 = 32 / ! + -g- ( v - !) + - 53 - v l v 

- I ) 2 

(182) 


and the 1 /rev critical region boundary on the y - u plane becomes 


= u 2 [- | (v 2 - 1) 3/2 ± \ /l + ^ (v 2 - 1 ) ♦ ^ v 2 (v 2 - l ) 2 J 

(183) 

The general behavior of the roots in this 1/rev critical region is 
correct but there are order y 2 errors in both the frequency and damping com- 
pared to the exact solution. C 2 , the width of the critical region, is exactly 
correct. There is an order u 2 change in the damping (eq. (179)); and the 
offset of the critical region, Cj, is not correct, leading to an order y 2 
change in the frequency. These order y 2 errors correspond to those in the 
expression for the roots away from the critical region (eq. (169)). For 
v * 1 the damping change goes to zero, and for v = 1 and Kp * 0 the offset 
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Ci reduces to zero for both the correct solution and for this approximation; 
in general, however, there are small errors in the 1/rev critical region 
behavior as there were for the roots away from the critical region. 

The coning mode (at 1/rev nonrotating frequency) and the high frequency 
mode (at 2/rev) do not participate at all in the 1/rev critical region. These 
roots are given still by their respective expressions for away from the 
critical regions (eqs. (164) and (169)). 


Summary 

The flapping dynamics of a rotor with three independent blades are 
described by a total of six eigenvalues. In the rotating frame, there are 
three independent equations; then the eigenvalues occur as two triple poles. 

In the nonrotating frame, the degrees of freedom and equations of motion are 
coupled. There are two eigenvalues at the rotating value for the coning mode, 
and four eigenvalues at the rotating values ± 1/rev for the high and low fre- 
quency rotor modes. Thus, if Xp is the rotating eigenvalue, the nonrotating 
roots are X^r * Xp, Xr ± i (for N = 3) . The total number of eigenvalues in 
the nonrotating frame is still six for the three-bladed rotor. The constant 
coefficient approximation in the nonrotating frame gives results for all six 
eigenvalues at once (in general, for all 2N roots of the rotor flap motion). 
This approximation gives the following behavior for the roots. 

Away from critical regions (i.e., the hover rotating frequency not near 
Vrev or 1/rev) the two coning mode roots are given by equation (164) , which 
is exactly the result obtained including the periodic coefficients. The four 
roots of the low frequency and high frequency modes are given by equa- 
tion (169), which has order y 2 errors in both the frequency and damping. 

The effects of forward flight away from the critical region are, in general, 
small, however, so these errors are not too significant. 

The four coning and low frequency mode roots encounter the Vrev critical 
region (at ±h/rev nonrotating frequency) , with exactly the same behavior as 
the correct solution (to order y at least) . The two roots of the high fre- 
quency mode (at 3/2/rev) do not participate in this behavior, which means an 
order y error. 

The four coning and high frequency mode roots do not encounter the 1/rev 
critical region, which means order y 2 errors for these roots. The two roots 
of the low frequency mode, at the real axis in the nonrotating frame do show 
the critical region behavior. There are, however, in general order y 2 
errors in both the frequency and damping compared to the correct solution, 
corresponding to the behavior of the roots of this mode away from the critical 
regions. The magnitude of these errors is discussed below in terms of the 
•, - y plane. 

Consider the y - y plane for v = 1 and Kp = 0 (fig. 2). The constant 
coefficient approximation has an identical 1/rev region boundary but only for 
the low frequency rotor mode, that is, two out of the six roots. It has the 
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same Vrev region, but only for the coning and low frequency modes (four out 
of six) . The real root boundary for the coning mode is the same but the 
boundary for the high and low frequency modes shows only half the influence of 
y. The latter boundary is, )n fact, only the point of closest approach of 
these two modes at ±l/rev nonrotating frequency, not an intersection at all. 

For v * 1.1 and Kp = 0 (fig. 4), the y - y plane of the constant 
coefficient approximation has a 1/rev region boundary, but only for the low 
frequency mode again. The region has the same width for a given y as the 
correct solution, but is offset higher. The error in the boundary is roughly 
Ay = 4y 2 , or about Ay = 1 at y = 0.5 (compared to y 0 about seven for the 
1/rev region with v = 1.1, as shown in fig. 4). As a result, rotors with y 
such that the hover frequency is above 1/rev will have a higher U C omer by 
this approximation. For example, with case (b) above, v = 1.1 and y = 6, the 
constant coefficient approximation gives ycomer = 0.326, compared to the 
correct solution of ycorner = 0.286. This error in the 1/rev region 
boundary - which corresponds to the frequency error away from the critical 
region, and so is a measure of the magnitude of that also - is not negligible, 
but is small in terms of the Ay shift of the boundary, or even in terms of 
the ycorner , which is more sensitive to the boundary shifts. The comparisons 
of the Vrev region and real root boundaries for this case follow exactly as 
for v = 1, discussed above. 

The constant coefficient approximation to the rotor flap dynamics in 
forward flight produces differential equations that do not actually describe 
the real rotor any more, and must always give some erroneous results. As a 
representation of the actual rotor dynamics, however, the constant coefficient 
approximation (in the nonrotating frame) is actually remarkably good. The 
influence of forward flight on the roots is given rather well by this approxi- 
mation, the primary error being that the roots of the high frequency rotor 
modes do not encounter the critical regions. For this case of N * 3, the 
Vrev region is seen by the four roots at ±Vrev, but not b/ the two roots at 
±3/2/rev; the 1/rev region is seen by the two roots at the real axis, but not 
by the four roots at ±l/rev and ±2/rev. This behavior is fundamental to the 
constant coefficient approximation; the roots from that approximation must 
always V complex conjugates, so the critical region behavior can only be 
exhibited by two roots on the real axis (as the 1/rev region here) ; or at a 
multiple of Vrev when there are four roots to participate, two at the posi- 
tive frequency and two at the negative frequency (as the Vrev region here). 
With these restrictions, the constant coefficient approximation picks up the 
critical region behavior of the periodic system whenever possible. The con- 
stant coefficient approximation in the nonrotating frame is better than in the 
rotating frame because the transformation to the nonrotating frame shifts the 
frequency of the roots to allow such occurrences. At fie least, these results 
suggest that the constant coefficient approximation is adequate for the low 
frequency dynamics of a rotor as for helicopter stability and control 
investigations. 
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Four-bladed Rotor 


Comparing the constant coefficient approximation of equations (156) and 
(157) for the N = 3 and N * 4 cases, it follows that the four-bladed rotor 
only adds the 62 degree of freedom and equation which are completely 
decoupled for the constant coefficient approximation. So, to equation (161) 
must be added the equation 


62 + } 62 + [v 2 + (1 ♦ U 2 ) J K p ]b 2 = 0 (184) 

for the four-bladed rotor. The roots for this mode are 

* * - je ± i/v 2 ♦ (l ♦ v 2 ) J k p - (X) 2 (185) 

This is the same as the constant coefficient approximation in the rotating 
frame which is not very good. All the critical regions are missed and there 
are even order p 2 errors away from the critical region. The solution for 
the other six roots is the same as for the N = 3 case. 


Five or More Blades 

As the number of blades increases, more and more modes are added to the 
nonrotating representation of the rotor. For example, N = 5 has degrees of 
freedom 8 0 , 8 lc , Bis. and 

62c = | cos 2 *m 

m 

&2S = | £l3 (m) Sin 2*n 
m 

with corresponding roots Xnr = Xr ± 2i. Consequently, more of the low fre- 
quency modes will be able to pick up the critical region behavior in the con- 
stant coefficient approximation. The highest frequency modes always do not 
encounter the critical regions, of course, so the constant coefficient approx! 
mation can never give completely correct behavior. This effect of increasing 
N parallels the effect on the differential equations. Increasing N tends 
to sweep the periodic coefficients from the lower frequency modes; there are 
always periodic coefficients present in the degrees of freedom and equations 
of the high frequency modes, however. 
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TRANSFER FUNCTIONS 


The experimental determination of the eigenvalues, especially the small 
increments (order y or even y 2 ) found here, in the harsh aerodynamic envi- 
ronment of the helicopter in forward flight, is a difficult task to accomplish 
with accuracy. A more direct measurement - hence, fundamentally more accu- 
rate - is the transfer function, that is, the response of the blade flap 
motion due to Sinusoidal excitation. 

Consider the transfer function of the flap response to blade pitch 
control, particularly the influence of the periodic coefficients in forward 
flight, which has led to the special behavior of the eigenvalues. The equa- 
tion of motion for an independent blade excited by pitch control inputs only, 
in the rotating frame, is 

6 + (J + J y sin 

+ I v 2 + ^ y cos i|i+ | u 2 sin 2\p + Kpf^- (1 + y 2 ) + ^ y sin y 2 cos 2i|/ll 6 


= ^ (1 + y 2 ) + j y sin iji - y 2 cos 0 (186) 

The transfer function is defined as the response to sinusoidal input 
6 = Fei^, where 0 is a complex constant. Taking only the real parts of 
both the input and the output is implied. With a constant coefficient differ- 
ential equation, the output will also be a sinusoid at frequency id, that is, 

6 = Fe^ 10 ^. The differential equation relates the output F to the input F 
by a single complex function H(id), which is the transfer function: 

1 = H(<d) (187) 

Q 

With a periodic coefficient differential equation, the response to 6 at 
frequency w is not just 6 at frequency id. The sinusoidal input at u> is 
multiplied by the coefficients which have terms that are also sinusoidal now 
with frequencies 1/rev, 2/rev, etc. (fl, 2Q, . . .). The product of two sinus- 
oidal functions is a sum of sinusoids at the sum and difference of the fre- 
quencies. It follows then that for a periodic coefficient system an input at 
frequency id leads to an output with terms at frequencies <d, u> ±l/rev, 
id ±2/rev, etc. The output is then a sum of the form 
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( 188 ) 


4-. . H,,!,). 11 " 11 ' * H , (u )e i( “- 1) * 

0 

+ H^uOe 1 ^* + H. 2 (u))e i(<1) - 2) * + . . . 


where H Q , H+j, H+ 2 , . . . are all the transfer functions ofc the system. 

Notice this may be written as 0/F = H(w, where H 0 , H ±1 , H+ 2 , etc. 

are the harmonics of a Fourier series representation of a function H(u>, i|>) , 
which is periodic in i|» (H(u), i|i) is complex, however, so H +n and H_ n are not 
conjugates). Equation (188) is a general lesult for periodic coefficient 
differential equations: the dynamic behavior is described by not a single 
transfer function, but rather a series of transfer functions. There is the 
direct response Hq, the output at the same frequency as the input; and there 
are also now sideband responses H+ n , output at the input frequency ±n/rev. 
Since for the hover limit y = 0 the present differential equation reduces to 
constant coefficients, all the transfer functions but H 0 must be zero then. 
It will be found, in fact, that H+ n are order u n . 


t / 
( 


Analysis 

Assume the input is 0 = ee 1 ^. It is possible to simply substitute into 
equation (186) the expansion for 8 as in equation (188), collect like har- 
monics, and thus solve for the transfer functions. Following the rest of the 
present investigation, however, an expansion in y will be used which makes 
the analysis more orderly. Expand the output as a series in y: 

8 = 6o + ligj + y 2 8 2 + . . . . Since it is known that the output has only time 

behavior like ei( w±n )'J , , there are no other time scales but i|i in this prob- 
lem. This means there is no critical region behavior in this problem that is 
a feature of the eigenvalues only. The critical region behavior is replaced 
here by the sideband transfer functions which will be shown to carry equiva- 
lent information about the system dynamics. It is therefore also not 
necessary to expand y as a m y. 

To order 1 the differential equation is 


Bo * o + + g - Kp^Bo = ^ 9 

So 8 1 gives Bq s Bo e ^ a) ^ with 
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where 


A(u) = -w 2 + |- iu) + v 2 + Kp | 


(191) 


The order 1 solution is the hover limit as usual; here it is the hover trans- 
fer function for the response of blade flapping to pitch input. The transfer 
function has the usual form for a second order system response in this case 
a heavily damped system with natural frequency (D n = v < (y/ 8 )Kp. Notice 
that A(iA) is the hover characteristic equation. 

The order y differential equation is 


gj + J Bi + (v 2 ♦ Kp J) 81 = J sin t |>6 - [j sin «J >& 0 ♦ (£ cos <|> ♦ * F \ sin <J»)b 0 ] 

= ^ l- 2 ie - (w + 1 - 2ikp)6oJe^ w 

. £ 1216 - (u,- 1 - ZiKpiBole 1 ^- 1 ^ 0«) 

Hence, the solution is 6 , ■ i,/'** 0 * * «l-‘ i( “‘ 1) * “ ith the tr “ sfer 
functions 


— " A (to ± l)A(u>) 


(193) 


The order y response is then just the u ± 1/rev sideband functions. The 
forcing terms, the right-hand side of equation (192), are the same terms which 
gave the order y, Vrev critical region for the eigenvalues. That is, near 
Vrev these terms contributed to the order y secular equation. Hence, the 
order y sideband transfer functions at w ± 1/rev constitute behavior of the 
system equivalent to the Vrev critical region of the eigenvalues. 


The order y 2 differential equation is 
h * } ♦ (v 2 + J Kp)B 2 = J (1 - cos 2H06 

- ^ sin 1 J/B 1 ♦ cos ♦ Kp j sin ^Bi 
♦ ^ sin 2\|) + Kp -g- - Kp ^ cos 2i(^BoJ (194) 
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Substituting for 0, Bg, and Bj, the solution is 


B 2 , B.e 1 "* . B 2 ,e 1( “* 2) * . Bz.e 4 


(195) 


with 


where 


«“> | - i - xrj 
8 


(ji 

A(<u) 


. y 2 (» - 2iKp)4i[v 2 - w 2 + ^(g + 2ia >)] 

+ (n) AfwHiCu) 


(u))Ai (to) 

Y ^2 (a) - 2iKp)4i(ui 2 - v 2 ) 
Ai((u) 


(12) 


(196) 


2 

Ai((o) = A(o) + 1 )A(j - 1) = [A (uj) - l] 2 ♦ ^2iw ♦ 


(197) 


Notice that A x (iX) is just the characteristic equation for the Bjc/^is 
hover roots. The direct response to order y 2 is then 


do 


W 2 B 2 


I 

8 

. 


(1 ♦ w 2 ) 


(12) 


2 (4iu) ♦ 8Kp) (to 2 - 


vfij 


v iw ♦ v 2 ♦ (1 + y 2 ) Kp 


* “ 2 


2 (4i(u ♦ 8Kp) |u 2 - v 2 ' '15 (g + 2iw)J 


‘1 


(198) 


So the order y 2 response contributes to the direct transfer function nuch as 
it does to the eigenvalues away from the critical regions. The order y 2 
correction for the mean of Mg is recognizable as usual, and there are other 
order y 2 effects due to the periodic coefficients. The u> i 2/rev sideband 
tr-r.jfer functions are 
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( 199 ) 


B2± g 1 

J * A(w i 2) 


Y , 2 (l6) (±i + Kp) - / Y N 2 ( “ 1 2 ' 2iK P } [' I 4 2i(( ° 2 ' v2) ] 
1A + A(u>) + V12/ A(u)A(w ± 1) 


Summary 

The transfer functions describing the response of the blade flap motion 
to pitch excitation are 


F = 


H±l ^ 


H 


±2 


8 0 + p 2 6 2 ' 


0 


2 ^2 + 
V -=r- 

0 


( 200 ) 


By measuring the response to a sinusoidal input, it is then possible to verify 
the equation representing the flap dynamics. Specifically, it is possible to 
verify the periodic coefficients which produce the critical region behavior of 
the eigenvalues. There is not, however, any critical region behavior for the 
transfer functions, that is, specific ranges of the parameters where there is 
more critical behavior of the response. That is replaced by the r* deband 
transfer functions. Hence, the transfer function measurement does not provide 
experimental demonstration of the critical region behavior, although it does 
equivalent informat’^n, verifying the equation which produces the critical 
region for the eigenvalues. 


A J'oint About Experimental Technique 

Since the response of a periodic system to input at frequency m is 
output at a) ± n/rev, for all integers n (although the magnitude of the 
response decreases with n) , it follows that the random excitation technique 
for measuring the transfer function of a constant coefficient system is not 
directly applicable to periodic coefficient systems. A random input has all 
frequencies at once, hence, the output is also composed of all frequencies. 

For constant coefficient systems the knowledge that the response at u came 
only from the input at m can be used to determine the transfer function when 
the input is random - the entire function at once, in fact. With a periodic 
coefficient system, however, the output at m has contributions from inputs 
at u ♦ n/rev for all n, so not enough information is available to determine 
the transfer function. 
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An input at a single frequency to must be used then and the output at 
all harmonics to ± n/rev measured; this ta c k is performed for the range of u 
required. This technique is slower than with the random excitation; in addi- 
tion, the noise filtering feature of the random input is lost, that is, 
accuracy of the measurements is lost. 

The problem is that there is not a single transfer function to measure 
but rather several. For example, the rotor flap response up to p = 0.5 or so 
would require the measurement of H 0 , H+i, and H+ 2 at least. It is probably 
possible to extend the random excitation techniques to periodic systems. It 
would be necessary, however, to use a number of inputs with independent 
spectra so that enough independent measurements are made at each frequency to 
determine the required number of transfer functions. In any case, the peri- 
odic system will require more experimental effort than the constant 
coefficient system. 


Nonrotating Response 

For the dynamics of the helicopter as a whole, it is the response of the 
rotor in the nonrotating frame that is of interest, that is, the response of 
the 8 lc» 8 ls> etc. degrees of freedom. Similarly, the inputs usually 

availablt for the rotor are collective and cyclic pitch control: 0q, 9ic 

^i s . He; ;, the transfer functions in the nonrotating frame may be examined, 
ui r.«.ditxon to the rotating response gi/en above. The analysis proceeds much 
Mve. It is straightforward since there are no critical regions to be 
concerned with or any other singular behavior; it is convenient to use vectors 
of the input and output variables and matrices of the transfer functions ($ 
due to ?) . The analysis and results are not very illuminating, however, 
because of their complexity. 


DISCUSSION OF PREVIOUS WORK 


Horvay (ref. 1) considered the flap dynamics of a rotor with v = 1, 
K p = 0, and no reverse flow eff~~ts. By the substitution 


6 = ye 


X 

16 




u cos 


* 


he transformed the flap equation to the standard form of Hill's equation: 


y + f(4»)y = 0 


where f is a periodic function. This equation he solved by the infinite 
determinant method of classical theory. The solution produced numerical 
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results presented as constant ReX lines on the y - u plane. The 
boundaries found were essentially those given in figure 7 although he did not 
investigate small enough y to find the 1/rev region. Horvay used the nota- 
tion n = y/8, since this is the coefficient of the flap damping in hover. 
Then, in hover and outside the critical regions, ReX = -y/16 = -n/2. Hence, 
he introduced the "apparent damping coefficient" n a pp = -2ReX. Then 
napp/n = -ReX/(y/16) takes the values 1 outside the critical regions, 1 to 0 
in the stable portions of the critical regions, and less than 0 in the 
unstable portions of the critical regions. The use of the notation n a pp/n 
and the parameter n a pp for ReX is found in much of the literature on rotor 
flap stability. Horvay found approximate solutions at high y (to y = 8) for 
several equations which are not true descriptions of the flap dynamics but do 
illustrate some of the high y behavior of the real rotor. 

Horvay and Yuan (ref. 2) considered blade flap stability for the cas j of 
v = 1 and no reverse flow but including Kp = 0. The roots were found '.rom 
the transient solution for 8, using the results of Floquet theory (as out- 
lined in Appendix A here). Much of the literature uses this approach, the 
real differences being the method used tc find the transient solution for 8. 
Horvay and Yuan used the ripple method to integrate the differential equation. 
They presented the numerical results on the y - y plane for v = 1 and 
Kp = 0, -/3/1Z, and /3/4 (these values of Kp were chosen because they place 
exactly at y/8 = /3 the h/rev region, the real root, and the 1/rev egion 
boundaries, respectively, when y = 0) . A large azimuth interval was used in 
the ripple method but the results compare fairly well with the solution of 
Horvay (ref. 1), at least at small y. 

Parkus (ref. 3) considered the flap stability for v = 1, no reverse flow 
but Kp t 0 (with the notation k = -Kp) . He found the roots from the Floquet 
result using an expansion of 8 as a series in y to find the transient 
solution. The expansion is the same as is used here away from the critical 
regions since the multiple time scales and expansion of y are not needed 
there. Parkus left his result in the form of a quadratic equation ~~>r 

a - e 2lt (y/ 16) ] . since it is only an order y 2 solution, howevn’ it is 
possible to solve explicitly for X. The result is exactly the same as 
obtained here for the roots away from the critical regions (eq. (24), with 
v = 1). Parkus did not recognize, however, that this solution is not valid 
for ImX near 1/rev or Vrev. The assumption that the 8 n in the expansion 
of 8 are ail the same order is violated near the critical regions. Parkus 
calculated the roots for varying y, for y = 12 and 14; the behavior of the 
roots looks like that in critical regions but really the expression for the 
roots away from the critical regions is breaking down as the roots approach 
the Vrev region. When ImXo is near the order y 2 term in the radicand 
of X (eq. (24) here) is large, so for large enough y the radicand is 
negative. The result is two real loots and an order y 2 effect which is far 
from he correct Vrev region behavior. 

Gessow and Grim (ref. 19) considered the flap stability at high y 
numerically integrating the equations of motion to find the transient behavior 
of 8. Reverse flow was included, as well as the effects of stall, compressi- 
bility, and large angles. They found the flap motion to be stable at y = 1 
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for y = 2.4 to 14.9, at u = 2.2 for y = 2.4 to 7.1, and at y = 3 for 
Y = 0.6 to 2.4. The motion was unstable at y = 3 for y = 14.9. 

Shulman (ref. 4) considered the flap stability for v - 1, Kp = 0, and no 
reverse flow. He added, however, a second degree of freedom, the first 
elastic flapwise bending mode of the articulated blade. His principal solu- 
tion method was simply to numerically integrate the equation of motion. He 
found a significant effect of the second degree of freedom for high advance 
ratio, y = 1.0 say, but small influence for y up to 0.5. The equations were 
solved for y = 1 and 1.2, but with the neglect of reverse flow the true flap 
dynamics are no longer represented at that y. His conclusion that the flap 
instability occurs at y = 1.5 to 1.8 is not correct then; but the importance 
of the second degree of freedom at high y is probably qualitatively correct. 

Perisho (ref. 5) considered the stability of the flap and pitch dynamics 

of the blade at high y, for an articulated rotor, Kp / 0, and including 

reverse flow aerodynamics. He solved the equations by numerical integration 
for the case v = 1.05, y = 5, and a pitch natural frequency of 10/rev. For 
just the rigid flap degree of freedom, he found the instability boundary at 
y = 2.2 with Kp = 0, and at y = 2.45 with Kp = 2, in the 1/rev region for 
both cases. Adding the torsion degree of freedom reduced the y for 
instability to y = 1.8 for Kp = 0; and to y = 1.65 for Kp = 2, in the 
%/rev region now. Adding flapwise bending, that is, three degrees of freedom, 
reduced the speed for instability still farther to y = 1.43 for Kp = 2. The 
reduction in the stability at high y due to the torsion degree of freedom 
appeared to be largely a torsion divergence in the reverse flow region. 

Shutler and Jones (ref. 6) considered the flap stability for v = 1 and 

Kp = 0, with no reverse flow. They solved for the eigenvalues and for 3 by 

a perturbation solution based on the Floquet theory result that the solution 
of a periodic coefficient differential equation may be written in the form 


8 = Cie^^ui^) + C2e^ 2 \i2(40 

where the eigenvectors uj and U 2 are periodic (Appendix A). The substitu- 
tion 8 = e^$u is made into the flap equation and then X and u expanded as 
series in y. The requirement that all the functions in the expansion of 
u(tJ0 must be periodic is equivalent to the requirement in the method of multi- 
ple time scales that the successive functions grow no faster than the earlier 
ones (the order 1 term of u(tJ0 is the hover limit, so it is constant; then 
requiring successive terms in the expansion grow no faster means they can at 
most be periodic); the result has the same secular terms as obtained here. 

This is basically the method used by Shutler and Jones, although the details 
differed considerably. To treat the critical regions it is necessary to 
quantify the requirement that the frequency be near a multiple of Vrev. The 
present investigation used an expansion of y to do this, while Shucler and 
Jones essentially expanded ImX itself (their eq. (25); with the present 
notation for X). So contrary to their statement, this expansion (eq. (25)) 
does have physical significance, namely, what "close" to a critical region 
means. Similarly, their parameter o also has physical significance: it 
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gives contours of constant ReX in the critical regions; for example, the 
critical region boundary is such a contour. They found the eigenvalues and 6 
for the Cc es of ^he %/rev region, the 1/rev region, and away from the criti- 
cal regior to order y 2 everywhere, including in the %/rev region. Their 
results m- ' be shown to be equivalent to those obtained here, for v = 1 and 
Kp = 0; except that here the %/rev region is only carried to order y. Their 
result for the %/rev region boundary to order y 2 has been given above 
(eq. (55a)) and illustrated in figure 7. They did not, however, take 
advantage of tho fact they had a perturbation solution to put together 
explicit analytic expressions for the eigenvalues. 

Lowis (ref. 7) considered the flap stability for v = 1 and Kp = 0, at 
high y particularly, fi.'st without and then with reverse flow effects. He 
solved for the roots by the results of Floquet theory, after numerically inte- 
grating the equation to find 8. Without reverse flow he found an instability 
at about y = 1.42 in the 1/rev region. With reverse flow, the instability 
occurred at about y = 2.3, still in the 1/rev region. The representation of 
the reverse flow aerodynamics used was only approximately correct, however: he 
did not account for the azimuth range where the blade is partly in normal flow 
and partly in reverse flow. This model has the correct limit for very large 
y but it is not good at ali below y = 1. The behavior of the solution 
around y = 2, that is, including the flap instability, is probably correct 
but neither of the solutions (without and with reverse flow) he obtained is 
good in the range y = 0.5 to y = 1.5 or so. 

Wilde, Bramwell, and Summerscales (ref. 20) considered flap stability at 
high y, including the effects of reverse flow, and presenting results mainly 
for v = 1 and Kr. = 0. They considered a teetering rotor also. The solution 
was obtained by use of an analog computer. They found a flapping instability 
at about u = 2.25 for y = 6. The teetering rotor was stable to y = 5 at 
least . 

Sissingh (ref. 8) considered the flap stability at high y, including 
v> 1, Kp / 0, and reverse flow (v/ith the notation P = v, Cj = -Kp; and only 
Kp = 0 was used for the results). The aerodynamic coefficients of reverse 
flow were derived and discussed. The solution of the equation of motion was 
obtained by use of an analog computer. The results were presented as stabil- 
ity boundaries for several y, on the y - v plane, with the emphasis on 
quite large v, typical of a slowed or stopping rotor. For v = 1, he found 
y = 2.2 for instability at y = 8 (higher for other values of y) . 

Stammers (ref. 9) considered the influence of forward flight on the 
flutter of the helicopter rotor blade. He found a perturbation solution by 
methods similar to those of Shutler and Jones (ref. 6) r< ther than the methods 
used in the present work. 

Hall (ref. 10) investigated the dynamics of the flap/lag/torsion motion 
of the blade. He used Floquet theory to obtain the roots from a numerical 
integration of the equations over one period. He discussed Floquet theory and 
it« results for a multi-degree of freedom system. He also presented results 
lor a single degree of freedom, that is, just rigid flap. 
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Sissingh and Kuczynski (ref. 11) extended the earlier work of Sissingh 
(ref. 8) to include the torsion degree of freedom; the emphasis again was on 
the results for very high y. They showed a quite significant effect of 
torsion on the blade stability reducing the y for instability considerably. 
For example, with y = 6 the y for instability was reduced to about 1.7 for 
a torsion frequency of 8/rev and to about y = 1.35 for a frequency of 5/rev. 
These boundaries were just a little below those for torsion only; the effect 
was largely due to torsional divergence in the reverse flow region where there 
is a negative aerodynamic spring and reduced aerodynamic damping on the pitch 
motion. Their results compared well with those of Perisho (ref. 5). 

Peters and Hohenemser (ref. 12) considered the flap stability including 
v > 1, Kp f 0, and reverse flow aerodynamics. They used a continuous repre- 
sentation of the reverse flow aerodynamic coefficients as obtained by Sissingh 
(ref. 8). They solved for the roots by integrating over one period and then 
using Floquet theory. Their solution compared well with the analog computer 
solution of reference 8. They presented the numerical results on the y - y 
plane for v = 1 and Kp = 0, 0.1, -0.1, out to y = 2.5. They found approxi- 
mately y = 2.3 for the flap instability in the 1/rev region (at y = 9, 
higher for other values of y) ; Kp > 0 increased the y for instability. 

There were possibly some numerical problems with the 1/rev boundary at low y 
(below 0.5), since the boundaries appear more like order y than the correct 
order y 2 behavior (their fig. 3). In addition, the y = 0 point of the 
1/rev region really must shift for Kp j- 0 and not remain always at y = 0 as 
shown (their figs. 4 and 5). 

Hohenemser and Yin (ref. 13) investigated the flap dynamics using the 
nonrotating degrees <>? freedom and equations including the periodic coeffi- 
cients in forward flight. They gave the equations for the case of N = 4. 

They solved for the riots by the methods of reference 12. They discussed the 

constant coefficient approximation to the nonrotating equations in forward 
flight but only in the context of tip path plane tilt feedback control, how- 
ever (at y = 0.4 with N = 3, ana y = 0.8 with N = 4) They found the con- 
stant coefficient approximation was not bad at low gain, especially for the 
low frequency modes but for high gain it could be unconservative. 

Johnson (ref. 14) considered rotor flap stability including v >. 1, 

Kp t 0, and reverse flow. He obtained perturbation solutions for the cases of 
small and large y and small and large y. The small y solution formed the 

basis for the analysis of the individual blade case in the present work. This 

report was summarized in reference 15, and reference 17 presented a synoptic of 
the small y results. 

Tong (ref. 16) considered the blade flap/lag dynamics. He obtained a 
solution by perturbation techniques to handle the nonlinear features of this 
problem (i.e., limit cycle instabilities), as well as the influence of forward 
f 1 ' ght . 

Biggers (ref. 18) considered the flap dynamics for Kp = 0 and no 
reverse flow but including v >. 1. He constructed y root loci and y - y 
planes for several cases using a numerical calculation of the exact 
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eigenvalues by means of Floquet theory. He presented the nonrotating 
equations of motion including the periodic coefficients in forward flight for 
M = 3 and N = 4. Then, he considered the constant coefficient approximation 
in the nonrotating frame, solving the characteristic equation for the roots of 
the system numerically and comparing this constant coefficient approximation 
with the exact results (Floquet theory) . He examined several cases out to 
y = 0.5. The results of this comparison of the exact and constant coefficient 
eigenvalues (both obtained by numerical techniques) are in agreement with the 
results of the comparison given here (of corresponding perturbation solutions). 


CONCLUSIONS 


This report has considered the influence of forward flight on the 
flapping stability of several helicopter rotor configurations. The eigen- 
values of the motion have been obtained by a perturbation technique which 
gives analytic expressions for the roots. Comparison between numerical solu- 
tions for the exact roots and the present perturbation solutions indicates 
that the latter are quantitatively accurate to about y = 0.5. (An exception 
is near the Vrev region, where the perturbation solution was carried to only 
order y. It should evidently be extended to order y 2 as was the rest of 
the solution.) In general, for this range of y, the flap motion retains the 
high aerodynamic damping of hover and so remains very stable in forward flight. 
There exist, however, critical regions due to the periodic coefficients which 
are encountered if the hover root frequency is too near a multiple of Vrev. 
The influence of the critical regions increases with y as the periodic 
coefficients increase. For the usual values of v, Kp, and y, the critical 
regions of interest are the Vrev and the 1/rev regions. In a critical region 
there is a plus and minus increment in the real part of the root from the 
hover value while the frequency remains fixed at a multiple of Vrev. Hence, 
there is a decrement in the stability of the system when a critical region is 
encountered. However, the damping change is only order y in the Vrev 
region, and order y 2 in the 1/rev region, so the stability decrease is small 
and the flap motion remains highly damped in forward flight. 

For y order 2 or so, that is, beyond the range or validity of the 
perturbation solutions obtained here, there can occur a sufficient stability 
degradation in a critical region (usually in the 1/rev region) so that a flap 
instability is encountered. 

These conclusions about the flap stability may generally be found in the 
existing literature on this problem. What the present work adds is explicit 
analytic expressions for the eigenvalues of the flap motion, including the 
periodic coefficient influence in forward flight. Also, the behavior of 
cantilever (v > 1) , teetering, and gimballed rotor configurations is examined 
in addition to that of an independent, articulated blade. 

The transfer function of the flap response to blade pitch has been 
considered as an alternative to the eigenvalues for describing the dynamic 
characteristics of the system. The transfer function indeed is found to 
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represent equivalent information about the dynamics, but with quite different 
behavior. The critical region behavior of the eigenvalues is replaced by 
sideband transfer functions (response at frequency <d ± nfl to input at 
frequency w) . 

The constant coefficient approximation to the flap equations of motion In 
the nonrotating frame was investigated and the eigenvalues compared with the 
solutions including the periodic coefficients. Such an approximation cannot 
be entirely correct, of course, but it is remarkably good, especially for the 
lower frequency nonrotating modes of the rotor. This implies that for certain 
problems - such as the low frequency dynamics of the rotor involving the heli- 
copter body motions - the constant coefficient approximation is an adequate 
representation of the system. The possibility of using the constant coeffi- 
cient approximation involves a considerable reduction in the effort required 
to analyze and understand the rotor dynamics. It is suggested, therefore, 
that the first step in an investigation involving helicopter rotor forward 
flight dynamics should be to check the validity of the constant coefficient 
approximation by comparing with an exact solution (obtained by numerical tech- 
niques probably) for the particular problem. The periodic coefficients (in 
the nonrotating frame) may not even be needed. If they are required, then the 
methods of perturbation theory are very useful in examining the fundamental 
behavior . 
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APPENDIX A 


PERIODIC SYSTEMS 


A. 
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The forward flight of the helicopter introduces periodic coefficients 
into the differential equations describing the flap motion due to the periodic 
variation of the free-stream velocity seen by the rotating blade. For large 
enough y the periodic aerodynamic forces radically influence the behavior of 
the root loci, and the analysis techniques required to find the eigenvalues. 

The root loci of a constant coefficient system characteristically exhibit 
behavior in which two roots start as complex conjugates, meet at the real axis, 
and then proceed in opposite directions along the real axis. The existence of 
periodic coefficients in the differential equations describing the motion 
generalizes this behavior so that it can occur at any frequency that is a 
multiple of %/rev, that is, at ImX = n/rev or n + %/rev where n is an 
integer, not just at 0/rev (the real axis) . The property cf the solution that 
allows this behavior is the fact the eigenvectors are themselves periodic 
(instead of constant as for a constant coefficient system; see the mathematics 
below) . The analysis which demonstrates that periodic systems show this 
behavior is called Floquet theory. 

Thus, the following behavior of root loci is characteristic of periodic 
systems (refer to fig. 10). If the parameter being varied, for example, the 
advance ratio y in the present problem, is such that at y = 0 the system is 
not periodic, then the roots start out as complex conjugates (point A on the 
loci in fig. 10). As y increases, the periodic forces increase, and the 
roots move toward n/rev (or n + %/rev) frequency, remaining complex conju- 
gates. At some critical y the loci reach ImX = n/rev (point B on fig. 10), 
and then for still larger y the frequency remains fixed at n/rev while the 
real pait of one root is decreased and that of the other is increased. The 
root being destabilized may cross into the rig 1 ’" half plane for some y 
(point C on fig. 10), indicating that the system has become unstable due to 
the influence of the periodic coefficients. 


A general system of differential equations with periodic coefficients may 
be reduced to a set of first-order equations, and may therefore be written 
(in matrix notation) as 

x = Ax (Al) 

where x is the state vector of the system and A(t) is a periodic matrix of 
coefficients: A(t + T) = A(t) . It may be shown that the solution to this 

differential equation can be obtained in the form 


x(t) = 

• .1 


A 


-i 
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7 * 


(A2) 





The Ai are the eigenvalues; the eigenvectors uj are periodic: 

u A (t + T) = Uj[(t); and the numbers q^(0) are constants obtained from the 
initial conditions. The theory that shows this is called Hoquet theory. The 
solution in this form is a direct extension of the normal solution for a con- 
stant coefficient differential equation which is characterized by constant 
eigenvectors (uj independent of time) . 

The eigenvalues A^ may be obtained by the following procedure. The 
equation 

<$> = A4> (A3) 

where <t> (t) is a matrix, is integrated over one period, from t = 0 to 
t = T, with the initial conditions 4>(0) = I (the unit matrix). Then, if A c ^ 
are eigenvalues of the matrix C = <f>(T) , the roots Aj are given by 
A c « or A = [In A c )/T. While the roots A c . (as eigenvalues of a real 

matrix C) must appear as real numbers or complex 1 con jugate pairs, the eigen- 
values A^ are under no such restriction, leading to the behavior of the root 
loci as described above. 

For a single degree of freedom, second-order system, let x R be the 
solution obtained from integrating the equation with the initial conditions 
x(0) = 1, x(0) = 0; and let xp be the solution with initial conditions 
x(0) = 0, x(0) = 1, Then, the roots A c are given by the quadratic equation: 

A c 2 - [x r (T) + x P (T)]A c + [x r (T)x p (T) - x R (T)x p (T) ] = 0 (A4) 
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APPENDIX B 


METHOD OF MULTIPLE TIME SCALES 


This appendix describes briefly the perturbation technique known as the 
method of multiple time scales. It should be used in parallel with the analy- 
sis of the flap equation, for many of the steps are more clearly expressed in 
the context of a specific example. More details of the method, and examples 
of its application may be found in reference 21 (and also ref. 14). 

Fundamental to perturbation methods is the existence of a small parameter, 
here the advance ratio p; and for the solution of periodic coefficient differ- 
ential equations, that the periodic terms be functions of p such that the 
equations reduce to constant coefficients when p = 0. For a study of the 
stability of a system, a solution is required that is uniformly valid over 
long time periods, so the long time behavior may be assessed. This leads to 
the use of the method of multiple time scales. Define a series of time scales 
i|/ n = y n ip . The time scales i|/ n are all assumed to be of the same order. Then 
for t|ij = piji the actual time i|» must be order p' 1 , that is, very large com- 
pared to the basic scale \1>q - The behavior of the solution over several 
time scales \p n will be investigated, each implying successively longer time 
behavior. The derivative with respect to time becomes then 

_9_ _ Jj_ _9_ 2 JL + 

9>|> 9||>0 ^ 9(|>i ^ 9^2 

So the assumption of the time scales is equivalent to an expansion of d/di(> 
as a series in p. 

Next the dependent variable is also expanded as a series in p: 


6 = <h» <l'2» • • •) + UBi (i|>o » ♦!» ...) + .. . 


The B n now depend on all the time scales i// n . The B n are all assumed to 
be of the same order, for all the long time scale behavior. This requirement 
is critical to obtaining the solution; it leads, for certain values of the 
free parameters, to critical regions in which there is a reduction in the 
stability of the system. In order to investigate the influence of the free 
parameters, they also are expanded as series in p. In this case, for the 
Lock number y have 


Y = Y 0 + UY 1 + W 2 Y 2 + • • • 


* 

i 



* 
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where the y n are all of the same order. For certain critical values of 
Y - that is, when y is order y or order y 2 from certain critical 
values - there will occur a stabilit) degradation due to the influence cf the 
periodic coefficients. 

So now 8, d/diji, and the free parameters are all expanded as series in y. 
These are substituted into the differential equation, and then all terms of 
like order in y are collected and separately set to zero. Thus, differ- 
ential equations are obtained for order 1, y, y 2 , ...» y n , .... With the 
expansion of d/di|>, that is, the use of multiple time scales, these equations 
are now partial differential equations for B n as functions of ij> n . The 
order y n equation has the following form: it may be written as a differ- 

ential equation for B n (t}i 0 ), forced by the lower order solutions. That is, on 
the right-hand side there are 8 n -l» 8n-2> ■ • • , Bo a nd their derivatives 
with respect to 4> n , ijj n _ j, . . <Pq. The set of partial differential equa- 
tions is solved progressively, beginning with the lowest order equation 
(order 1), then order y, etc. 

It is the characteristic of this expansion in y that the y n order 
equations, written as ordinary differential equations for Bn(’J’o) » all have 
the same homogeneous solution. Since the equation for BnC'J'o) forced by 
the lower order solutions 8 n -i» • • •> Bo, it follows then that when the 
known solutions for 8 n _j, . . . , Bo are substituted, the equation for 8 n 
will be forced by its own homogeneous solution. This would give rise to 
solutions for B n of the form 4>o times its homogeneous solution, that is, 
solutions for B n of order d»o^n- 1 » • • •» Then S n /6n-i» • • •» 6 n /Bo 

would be order ij'o, that is, would become arbitrarily large if <i>o is large 
enough. This violates the assumption that Bn is of the same order as 
B n -i» • • • , Bo- The only way such a solution may be avoided is if the coef- 
ficient of the homogeneous solution on the right-hand side of the y n order 
equation for BnC’i'o) is itself set to zero. This coefficient of the homogene- 
ous solution forcing the equation for 8 n (<J>o) is called the secular term. The 
secular term is set to zero so that the solution for B n will be uniformly 
valid for all time, that is, for large ij>o- 

Now the right-hand side of the 8 n equation itself involves derivatives 
of B n -i» • • - , Bo with respect to the time scales i|/ n , . . ., At this 

stage in the analysis, considering the y n order equation, some of the solu- 
tion for B is known already. What is not known is the behavior of Bn-i 
with respect to i|»i, \p 2 , . . . , the behavior of B n -2 with respect to 
i|> 2 , . . and so on down to Bo with resoect to i|» n . Hence, these combina- 
tions of the dependent and independent variables remain on the right-hand side 
of the B n (<|>c) equation when the known solution is substituted. Specifically, 
these combinations will be in the coefficients of the homogeneous solution, 
that is, in the secular term. Then, setting the secular term to zero results 
in a differential equation, which may be written as an equation for 
forced on the right-hand side by the lower order solutions. This equation 
will also be forced by its own homogeneous solution, and so its secular term 
must be set to zero in order to maintain the uniform validity of 8 n -i over 
the time scale. The result is an equation for B n _ 2 (^2)l this process 

is continued down to an equation for 3 q C4-n) - 


91 



So the order p n equation produces a set of differential equations for 
3 n 0 M» 8 n _j0h), S n -?^2)» • • •» Bo('i'n) • This set of equations is then 
solved, for the behavior of g n with respect to ipQ, B n -i with respect to 
i|>l, and so on down to 80 with respect to i|> n , thus completing the solution 
to order p n . The analysis then proceeds to the order y n+1 equation. For 
the stability of the system, it is the behavior of with respect to i|/ n 
that is of interest since that giv's the eigenvalue to order p n . The nesting 
behavior of the differential equations makes each successive order solution 
more involved but the basic procedure remains the same. 
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APPENDIX C 


SOLUTION OF THE SECULAR EQUATION 


The application of the method of multiple time scales to periodic 
coefficient differential equations results in secular equations of the form 


+ (a ♦ id) (5 ♦ (b ♦ ic)B = 0 


(Cl) 


where 6 is a complex quantity, and the constants a, b, c, and d are real. 
Writing 


D 2 = d 2 - (b 2 + c 2 ) = d 2 - |b + ic | 2 


(C2) 


D = /|D 2 | 


(C3) 


it may be verified that the solution of this equation is, for D 2 > 0: 
6 * e ai ^ | a [ d - D + i(b + ic) + A[d ♦ D ♦ i(b + ic) Je"* 1 ^ | 

where A is a complex constant; for D 2 = 0: 


(C4) 


6 * e'^CAtfd ♦ i(b * ic)]<|; ♦ i) ♦ B[d ♦ i(b ♦ ic)]) 


(C5) 


where A and B are real constants; and for D 2 < 0: 

’ (a [ d ♦ iD * i (b ♦ ic) ] e 0 * + B [d - iD ♦ i (b ♦ ic) ] e' 0 * J (C6) 


0 = e-^i 


where A and B are real constants. The limiting case b 
D = d, so the solution is 


c = 0 gives 


B ’ Ae 


- (a+idH 


(C7) 


where A is a :omplex constant. 

The regio.. of decreased stability - the critical region * is given by 
D 2 < 0. The boundary of the critical region is P 2 = 0. The behavior of the 
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solution of this equation when b or c i 0 _is that described for periodic 
systems; indeed it will be found that the $ turm comes from the periodic 
coefficients. 
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Figure 6.- Critical region and real root boundaries on the y - y plane for 
v = 1 and Kp * 0.1 (based on the perturbation solution). 
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Figure 8.- Critical region, real root, and high y instability boundaries on the y - y plane for 
v * 1 and Kp » 0, for y • 0 to 2.5 (a composite plot based on the results of refs. 1 to 18). 
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Figure 9.- Behavior of the y root locus of the three-bladed gimballed rotor, 
near the two real root boundary. The thr«e cases are discussed in the 
text. 








